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Abstract 


A  high  fidelity  simulation  of  the  coherent  signal  samples  produced  by  a  low-PRF  shipbome 
pulsed  Doppler  radar  when  observing  a  sea-skimming  anti-ship  missile  moving  through 
precipitation  above  the  sea  has  been  created.  The  modelled  radar  system  uses  a  phased  array 
antenna  and  has  the  ability  to  vary  the  parameters  of  the  transmitted  waveform  from  burst  to 
burst.  The  output  of  the  simulation  is  designed  for  combination,  with  the  output  of  clutter  and 
signal  processor  simulations  so  that  target  detection  in  clutter  can  be  examined.  A  complete 
description  of  the  physical  models  used  in  the  simulation  and  their  mathematical 
implementation  is  presented.  Example  output  data  were  used  to  derive  propagation  factor 
values,  which  were  then  compared  with  the  output  of  the  TERPEM  commercial  software 
package.  Received  power  values  were  also  derived  from  the  example  data  and  compared  with 
the  output  of  a  multifunction  radar  simulator  (ADAPT_MFR).  Good  agreement  between  the 
data  sets  was  achieved,  except  for  expected  differences  found  at  ranges  near  the  radar  horizon 
(“intermediate  region”).  Thus,  this  simulation  is  capable  of  producing  realistic  and  useful 
signals  for  target  detection  studies.  It  is  recommended  that  the  simulation  be  upgraded  by 
integrating  the  TERPEM  software.  This  will  allow  propagation  to  be  correctly  considered 
near  and  beyond  the  radar  horizon,  thus  extending  the  applicability  of  the  simulation. 

Resume 


Une  simulation  haute  fidelite  d’echantillonnage  de  signaux  coherent  produit  par  un  radar 
Doppler  naval  a  basses  frequences  de  repetition  de  pulses  observant  un  missile  anti-navire  qui 
se  deplace  au-dessus  de  la  mer  a  travers  des  precipitations  a  ete  creee.  Le  modele  de  systeme 
radar  utilise  une  antenne  a  reseaux  de  phase  et  peut  changer  les  parametres  de  la  forme  d’onde 
transmise  de  groupe  de  pulses  en  groupe  de  pulses.  Les  sorties  de  la  simulation  sont 
choisies  pour  se  combiner  avec  les  sorties  de  fouillis  et  les  simulations  de  traitement  de 
signaux  de  fatjon  que  la  detection  de  cible  dans  un  fouillis  puisse  etre  etudier.  Une  description 
complete  des  modeles  physiques  utiliser  dans  la  simulation  et  leur  representation 
mathematique  est  presentee.  Des  donnees  de  sorties  ont  ete  utilisees  pour  deriver  des  valeurs 
de  facteur  de  propagation,  lesquelles  ont  ete  comparees  avec  les  sorties  du  programme 
commercial  TERPEM.  Les  valeurs  de  pouvoir  requs  ont  ete  aussi  deliver  de  ces  donnees  et 
compare  avec  les  sorties  du  simulateur  de  radar  multifonction  (ADAPT_MFR).  Un  bon 
accord  entre  ces  groupe  de  donnees  a  ete  produit,  exceptee  pour  une  difference  prevu  au 
distance  proche  de  1’ horizon  radar  (region  intermediate).  Ainsi,  cette  simulation  est  capable 
de  produire  des  signaux  realiste  et  utile  pour  des  etudes  de  detection  de  cibles.  II  est 
recommande  que  la  simulation  soit  ameliorer  en  integrant  le  logiciel  TERPEM.  Cela 
permettra  a  la  propagation  d’etre  considerer  correctement  proche  et  au  dela  de  1’ horizon  radar, 
ainsi  elargissant  les  applicability  de  cette  simulation. 
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Executive  Summary 


Recent  emphasis  on  naval  operations  in  the  littoral  zone  and  the  potent  threat  of  sea-skimming 
anti-ship  missiles  has  increased  the  importance  of  understanding  the  capability  of  radar  to 
detect  small  targets  in  highly  cluttered  environments.  Simulation  is  a  useful  approach  for 
developing  an  initial  assessment  of  radar  detection  performance,  since  acquiring  real 
shipbome  radar  measurements  of  sea-skimming  missiles  in  littoral  clutter  is  a  major 
undertaking.  In  addition,  it  would  be  very  difficult  to  acquire,  experimentally,  a  data  set 
encompassing  the  full  range  of  clutter  conditions  that  shipbome  radar  must  operate  in. 

This  report  describes  a  high  fidelity  simulation  of  the  received  signals  produced  by  shipbome 
radar  when  observing  a  sea-skimming  anti-ship  missile  moving  through  precipitation  above 
the  sea.  The  simulation  was  designed  as  a  component  of  a  larger  simulation  aimed  at 
determining  the  detection  performance  of  radar  operating  in  clutter  conditions.  The  physical 
models  used  in  the  simulation  and  their  mathematical  implementation  is  described. 

This  report  also  presents  an  evaluation  of  the  target  simulation  by  comparing  quantities 
derived  from  example  output  data  with  similar  data  generated  by  a  commercial  software 
package,  and  by  a  lower  fidelity  multifunction  radar  simulation.  Good  agreement  was 
achieved  in  each  comparison,  except  where  the  calculations  were  extended  to  the  radar 
horizon.  It  is  recommended  that  this  simulation  be  integrated  with  the  commercial  software 
package  TERPEM,  so  that  its  applicability  can  be  extended  to  any  range  near  or  beyond  the 
radar  horizon  and  to  any  atmospheric  structure.  Nevertheless,  the  favourable  comparisons  in 
the  interference  region  demonstrate  that  this  simulation  is  accurate,  where  applicable,  and  that 
it  can  be  used  with  confidence  for  its  stated  purpose  of  assessing  radar  detection  performance 
in  clutter. 


Thomson,  A.  D.,  2002.  A  Target  Simulation  for  Studies  of  Radar  Detection  in  Clutter. 
DRDC  Ottawa  TR  2002-145.  Defence  R&D  Canada  -  Ottawa. 


DRDC  Ottawa  TR  2002-145 


iii 


Sommaire 


L’emphase  rdcente  sur  les  operations  navales  en  region  littorale  et  la  menace  existante  des 
na vires  anti-navire  a  augmente  1’ importance  de  comprendre  la  capacite  des  radars  pour 
detecter  des  cibles  petites  dans  des  fouillis  important.  La  simulation  est  une  approche  utile 
pour  developper  une  evaluation  initiate  de  la  performance  de  detection  radar,  puisque 
l’acquisition  de  donnde  rdelle  de  mesures  de  radar  navales  de  missile  anti-navire  volant  dans 
un  fouillis  littoral  est  un  effort  considerable.  De  plus.  D  serait  tr£s  difficile  d’obtenir 
experimentalement  un  ensemble  de  donnees  incluant  la  grande  variete  de  conditions  de 
fouillis  qu’un  radar  naval  peut  rencontrer. 

Ce  rapport  decrit  une  simulation  haute  fideiite  d’echantillonnage  de  signaux  re§us  par  un 
radar  naval  qui  observe  un  missile  anti-navire  se  ddpla^ant  a  travers  des  precipitations  au 
dessus  de  la  mer.  La  simulation  a  ete  con$ue  comme  une  composante  d’un  plus  gros 
simulateur  pouvant  determiner  la  performance  de  detection  de  radar  operant  dans  des 
conditions  de  fouillis.  Une  description  complete  des  modules  physiques  utiliser  dans  la 
simulation  et  leur  representation  mathemadque  est  decrite. 

Ce  rapport  presente  aussi  une  evaluation  de  la  simulation  en  comparant  des  quantites  ddrivees 
d’exemple  spedfique  avec  des  quantites  similaires  generer  par  un  logiciel  commercial,  et  un 
simulateur  de  radar  muldfonction  basse  fideiite.  Un  bon  accord  a  ete  produit  pour  chaque 
comparaison,  excepte  pour  les  calculs  au-dete  de  l’horizon  radar.  II  est  recommande  que  la 
simulation  soit  amdliorer  en  integrant  le  logiciel  TERPEM,  ce  qui  permettra  d’dtendre  ses 
applicabilites  k  des  distances  proche  ou  au  delS  de  1’ horizon  radar  et  aussi  pour  n’importe 
quelle  structure  atmospherique.  Cependant,  les  comparaisons  favorables  dans  la  region 
d’ interference  intermediate  demontrent  que  cette  simulation  est  precise,  quand  applicable,  et 
qu’elle  peut  etre  utiliser  avec  confiance  pour  revaluation  de  performance  de  detection  de  cible 
en  fouillis. 
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Figure  1.  The  general  scenario  modelled  by  the  target  simulation.  Blue  shading  depicts  the 
spherical  sea  surface.  Grey  shading  depicts  the  three  phases  of  precipitation  in  the 
atmosphere.  The  thin  black  line  depicts  the  missile  trajectory  at  a  constant  altitude  along 
a  constant  azimuth  with  respect  to  the  radar . . . 2 

Figure  2.  Overview  of  the  computational  approach  used  to  generate  the  target  signals.  Here 
the  dimension  of  the  abscissa  in  the  power  and  phase  spectra  plots  is  Doppler  velocity,  v, 
as  opposed  to  Doppler  shift  frequency . 3 

Figure  3.  Example  height  profile  of  reflectivity  factor,  Z,  where  dBZ  =  10.0  log(Z  in  mm6  m"3)7 

Figure  4.  The  range  weighting  function . . . 1 1 

Figure  5.  Antenna  geometry . 12 

Figure  6.  Azimuthal  scan  pattern.  The  depicted  3  dB  beamwidth  on  boresight  is  10°.  The 
azimuth  extent  of  the  scan  area  is  90°.  The  dashed  lines  show  the  boundaries  of  the  scan 


area  and  the  antenna  azimuth  boresight.  Grey  shaded  regions  depict  the  individual 
beams.  The  beam  for  which  radar  signals  are  calculated  is  shaded  yellow . 15 

Figure  7.  Geometry  of  the  missile  trajectory.  The  radar  antenna  is  located  at  the  origin  of  the 
y  and  z  axes.  The  thick  blue  curve  represents  the  surface  of  the  sea . 17 

Figure  8.  Example  of  burst  missile  positions  relative  to  the  radar  resolution  volumes . 19 

Figure  9.  Geometry  describing  target  positions  below  the  radar  horizon.  The  thick  blue  curve 
represents  the  sea  surface. . . 21 

Figure  10.  The  signal  paths  that  cause  interference.  The  thick  blue  curve  represents  the  sea 
surface . 22 

Figure  11.  Path  3  geometry . 24 


Figure  12.  Comparison  of  the  propagation  factors  calculated  by  the  target  simulation  and 
TERPEM.  The  modelled  scenario  is  represented  by  the  parameters  listed  in  Tables  1, 2, 
and  3.  The  Ts,  ft,  and  Nc  parameters  for  this  case  represent  burst  0  of  a  four  burst  dwell 
that  is  repeated  continuously  throughout  the  flight  time  of  the  missile.  The  blue  curve 

shows  Fj  vs.  range.  The  red  curve  shows  F'4  vs.  range.  The  green  curve  shows 

F4/Pc  vs.  range . 32 

Figure  13.  Comparison  of  the  propagation  factors  calculated  by  the  target  simulation  and 
TERPEM.  The  target  simulation  data  corresponds  to  the  burst  1  waveform  parameters. 
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The  blue  curve  shows  Fr4  vs.  range.  The  red  curve  shows  F'4  vs.  range.  The  green 
curve  shows  F4jPc  vs.  range . . . 34 

Figure  14.  Comparison  of  the  propagation  factors  calculated  by  the  target  simulation  and 
TERPEM.  The  target  simulation  data  corresponds  to  the  burst  2  waveform  parameters. 

The  blue  curve  shows  Fr4  vs.  range.  The  red  curve  shows  F'4  vs.  range.  The  green 

curve  shows  F4/Pc  vs.  range . 35 

Figure  15.  Comparison  of  the  propagation  factors  calculated  by  the  target  simulation  and 
TERPEM.  The  target  simulation  data  corresponds  to  the  burst  3  waveform  parameters. 

The  blue  curve  shows  Fr4  vs.  range.  The  red  curve  shows  F'4  vs.  range.  The  green 
curve  shows  F4/Pc  vs.  range . 36 

Figure  16.  Comparison  of  the  propagation  factors  calculated  by  the  target  simulation  and 
TERPEM.  The  parameters  of  the  modelled  scenario  are  the  same  as  that  corresponding 
to  Figure  15,  except  for  antenna  height,  which  was  set  to  5  m.  The  blue  curve  shows 

Fj  vs.  range.  The  red  curve  shows  F'4  vs.  range.  The  green  curve  shows  F4/Pc  vs. 
range . 37 

Figure  17.  Comparison  of  the  burst  0  target  powers  generated  by  the  target  simulation  and 
ADAPT_MFR  for  the  scenario  defined  in  Tables  1,  2,  and  3.  The  green  curve  shows  the 
standard  output  of  the  target  simulation.  The  red  curve  shows  the  data  generated  by  the 
target  simulation  without  range  weighting  or  statistical  fluctuations.  The  blue  curve 
shows  the  equivalent  target  powers  generated  by  ADPT_MFR . 37 

Figure  18.  Comparison  of  the  burst  1  target  powers  generated  by  the  target  simulation  and 
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the  standard  output  of  the  target  simulation.  The  red  curve  shows  the  data  generated  by 
the  target  simulation  without  range  weighting  or  statistical  fluctuations.  The  blue  curve 
shows  the  equivalent  target  powers  generated  by  ADPT_MFR . 38 
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the  standard  output  of  the  target  simulation.  The  red  curve  shows  the  data  generated  by 
the  target  simulation  without  range  weighting  or  statistical  fluctuations.  The  blue  curve 
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1  Introduction 


Recent  emphasis  on  naval  operations  in  the  littoral  zone  and  the  potent  threat  of  sea-skimming 
anti-ship  missiles  has  increased  the  importance  of  understanding  the  capability  of  radar  to 
detect  small  targets  in  highly  cluttered  environments.  Simulation  is  a  useful  approach  for 
developing  an  initial  assessment  of  radar  detection  performance,  since  acquiring  real 
shipbome  radar  measurements  of  sea-skimming  missiles  in  littoral  clutter  is  a  major 
undertaking.  In  addition,  it  would  be  very  difficult  to  acquire,  experimentally,  a  data  set 
encompassing  the  full  range  of  clutter  conditions  that  shipbome  radar  must  operate  in. 

This  report  describes  a  high  fidelity  simulation  of  the  coherent  signal  samples  produced  by 
shipbome  low-PRF  radar  when  observing  a  sea-skimming  anti-ship  missile  moving  through 
precipitation  above  the  sea.  The  simulation  was  designed  as  a  component  of  a  larger 
simulation  aimed  at  determining  the  detection  performance  of  radar  operating  in  clutter 
conditions.  The  physical  models  used  in  the  simulation  and  their  mathematical 
implementation  is  described.  These  models  are  the  same  as  or  compatible  with  the  models 
used  in  the  other  components  of  the  larger  detection  simulation.  In  some  cases  these  models 
represent  upgraded  versions  of  models  previously  described  for  the  weather  clutter  simulation 
component  [1].  In  such  cases  the  upgraded  models  have  also  been  applied  to  the  weather 
clutter  simulation.  Hence,  this  report  serves  as  a  description  of  upgrades  to  the  weather  clutter 
simulation.  This  report  also  presents  an  evaluation  of  the  target  simulation  by  comparing 
propagation  factor  values  derived  from  example  output  with  data  generated  by  the  commercial 
software  package  TERPEM.  In  addition,  the  received  power  data  generated  by  the  target 
simulation  are  compared  with  similar  data  generated  by  a  multifunction  radar  simulation 
(ADAPT_MFR,  Jones  et  al.  [2]). 
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2  General  Approach 


The  geometry  of  the  general  scenario  that  is  modelled  by  the  target  simulation  is  depicted  in 
Figure  1.  A  sea-skimming  missile  flies  at  constant  altitude  above  the  sea  directly  toward  the 
radar  system  attempting  to  detect  it.  The  azimuth  of  the  missile  approach  also  remains 
constant  with  respect  to  the  radar.  This  corresponds  to  a  worst-case  detection  scenario  in 
terms  of  the  radar  cross-section  of  the  missile.  The  simulation  could  be  expanded  to  include 
manoeuvring  targets  if  the  radar  cross-section  of  the  target  was  known  as  a  function  of  aspect 
angle.  Figure  1  also  shows  that  the  missile,  which  moves  at  a  constant  velocity,  has  a 
different  position  each  time  it  is  illuminated  by  a  radar  pulse. 

The  modelled  atmosphere  is  composed  of  gases,  clouds,  and  precipitation  that  will  attenuate 
the  signals  transmitted  by  the  horizon  scanning  radar.  Clutter  returns  and  radar  system  noise 
are  generated  by  the  other  simulations  that  makeup  the  overall  detection  simulation  (e.g.,  see 
Thomson  [3]).  Propagation  is  currently  considered  in  a  simplistic  manner.  Refraction  of  the 
radar  signal  is  considered  using  an  effective  earth  radius.  Interference  caused  by  reflection  of 
the  signal  from  the  sea  surface  is  calculated  geometrically.  Diffraction  by  the  sea  surface  is 
neglected,  which  places  important  limitations  on  the  applicability  of  this  simulation. 
Specifically,  this  approach  can  only  be  used  in  the  “interference  region”  (see  Blake  [4], 
p.  269-279).  When  the  modelled  target  is  within  the  “intermediate  region”,  near  the  radar 
horizon,  the  limited  propagation  calculations  will  be  incorrect.  This  is  not  considered  to  be  a 
major  problem  when  investigating  radar  detection  of  targets  within  strong  clutter  since  in 
most  cases  the  clutter  will  prevent  detection  of  the  target  in  the  intermediate  region  anyway. 
However,  care  should  be  used  when  analyzing  detection  results  to  ensure  that  target  detections 
within  the  intermediate  region  are  properly  considered. 
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Figure  1.  The  general  scenario  modelled  by  the  target  simulation.  Blue  shading  depicts  the  spherical 
sea  surface.  Grey  shading  depicts  the  three  phases  of  precipitation  in  the  atmosphere.  The  thin  black 
line  depicts  the  missile  trajectory  at  a  constant  altitude  along  a  constant  azimuth  with  respect  to  the 

radar. 
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The  main  outputs  of  the  simulation  are  the  discrete  in-phase,  \n ,  and  quadrature,  Q„ , 
components  representing  the  voltage  signal  received  by  the  radar  when  observing  the  target. 

In  and  Qn  values  are  calculated  for  the  range  bins  of  each  scan  in  which  the  missile  is 
illuminated  by  the  radar.  These  range  bins  will  only  be  found  in  the  radar  beam  that  covers 
the  constant  azimuth  of  the  missile  approach. 

Figure  2  presents  an  overview  of  the  mathematical  approach  used  to  calculate  the  target 
signals.  Although  the  simulation  produces  signals  for  many  range  bins  in  many  different 
scans,  the  figure  only  depicts  the  process  for  one  range  bin.  The  process  begins  with  a  missile 
model  and  a  statistical  model. 

The  missile  model  calculates  the  mean  in-phase,  I„  ,  and  quadrature,  Q„ ,  components  of  the 

target  signal.  These  calculations  consider  the  transmit  and  receive  characteristics  of  the  radar, 
attenuation  of  the  signal  by  the  precipitation  field,  the  target  characteristics,  the  sea  surface 
conditions,  and  multipath  interference  caused  by  a  reflective  sea  surface.  Details  of  these 
calculations  are  given  in  Section  5. 

The  modelled  missile  can  be  considered  to  have  several  different  scattering  features. 
Consequently,  as  the  orientation  of  the  missile  changes  with  respect  to  the  radar,  the 
interference  of  the  waves  backscattered  from  the  different  scattering  features  will  vary.  The 
statistical  model  creates  random  numbers  that  represent  the  time  varying  fluctuations  in  the 
radar  signal  that  result  when  the  interfering  waves  are  summed  at  the  antenna  port.  The  rate 
of  fluctuation  and  the  statistical  properties  of  the  random  numbers  output  by  the  statistical 
model  correspond  to  the  four  Swerling  models  (see  Skolnik  [5]).  The  particular  Swerling 


Figure  2.  Overview  of  the  computational  approach  used  to  generate  the  target  signals.  Here  the 
dimension  of  the  abscissa  in  the  power  and  phase  spectra  plots  is  Doppler  velocity,  v,  as  opposed  to 

Doppler  shift  frequency. 
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model  used  is  an  input  parameter  specified  by  the  user.  In  addition,  the  user  may  also  specify 
a  non-fluctuating  target  model  in  which  the  overall  I„  and  Q„  output  of  the  target  model  is 

taken  directly  from  the  missile  model,  i.e.,  I„  =  I„  and  Qn  =  Q„ .  Specifically,  the  statistical 

model  creates  vectors  of  random  numbers,  each  Nc  elements  long,  for  each  range  bin  where 
the  target  signal  is  non-zero.  Here  Nc  denotes  the  number  of  coherent  pulses  per  burst 

transmitted  by  the  radar.  The  two  sets  of  random  numbers,  f5=(f0s,f15,  ...,f,s,  ...,f^_j)  and 

f  *  =  (f  * , f* ,  ...,  f®, .. .,  f  *_! ),  represent  the  fluctuations  in  the  power  and  phase  spectra, 
respectively,  of  the  target  signal.  Details  of  the  statistical  model  are  given  in  Section  6. 

The  coefficients  of  the  mean  power  and  phase  spectra,  S,  and  O, ,  of  the  missile  signal 
corresponding  to  a  particular  range  bin  are  calculated  from  and  Q„  using  a  discrete 
Fourier  transform  (DFT),  i.e., 


5,=(Re{x,})!  +  (lm{x,})2 


<D,  =  tan"1 


Wx.iJ 


(1) 


where 

X  =  DFT{!  +  yQ}  j  =  (2) 

and  Re  and  Im  take  the  real  and  imaginary  part  of  their  respective  arguments.  The  signals  I 
and  Q  also  contain  Nc  elements  each,  i.e., 


Q  —  ( Qo  *  Qi » •  •  •  *  Qn » •  •  •  *  Qwc-i ) 


(3) 


The  statistical  fluctuations  are  incorporated  into  the  simulated  target  signal  by  filtering  f  5 
with  S  ,  i.e.,  the  power  spectrum  coefficients  are  given  by 

S,=S,  fts  (4) 


The  phase  spectrum  of  the  final  target  signal  is  given  by  the  mean  phase  spectrum,  <t> ,  for 
Swerling  cases  1  and  3,  and  is  given  by  the  phase  fluctuations,  f  ® ,  for  Swerling  cases  2  and 
4.  The  final  I„  and  Q„  signal  samples  are  calculated  by  inverse  DFT,  i.e., 
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1  = 


Q 


=  Im 


where 


X,  =  t[s~  exp  (j  0; )  for  Swerling  1  and  3  cases 

X;  =  t[s~ exp (;' /j® )  for  Swerling  2  and  4  cases 


(6) 
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3  Environment  Model 


The  modelled  environment  is  the  same  as  that  used  in  the  weather  clutter  simulation  that  also 
forms  a  part  of  the  larger  detection  simulation  referred  to  earlier  (e.g.,  see  Thomson  [3]).  A 
complete  description  of  the  modelled  environment  and  its  effects  on  the  radar  signal  are  given 
in  Thomson  [1].  Only  a  brief  summary  will  be  given  here. 

The  modelled  environment  consists  of  a  stratiform  precipitation  field  above  a  rough  spherical 
sea  surface.  The  precipitation  field  consists  of  horizontally  homogeneous  layers  of  rain, 
mixed  phase  particles,  and  snow.  A  uniform  cloud  field  also  exists  from  a  specified  lower 
level  to  the  top  of  the  snow  layer.  The  input  parameters  required  to  specify  this  environment 
are  listed  in  Table  1. 

The  quantity  Z(h)  shown  in  Table  1  denotes  the  height  profile  of  reflectivity  factor  where  h  is 
height  above  sea  level.  The  reflectivity  factor  is  related  to  the  backscatter  cross  section  of  the 
precipitation  (see  Doviak  and  Zmid  [6])  and  provides  a  measure  of  the  precipitation  intensity. 
For  example,  Z  is  often  related  to  the  rain  rate,  R,  through  a  power  law,  i.e.,  Z  =  aRh .  The 
input  Z  profile  is  extrapolated  horizontally  to  apply  to  any  location  within  the  viewing  space 
of  the  radar.  An  example  profile  is  shown  in  Figure  3.  In  this  example,  the  snow  region  is 
characterized  by  decreasing  Z  above  4.0  km  altitude.  The  region  of  larger  Z  between  3.4  and 

Table  1.  Parameters  describing  the  environment. 


SYMBOL 

EXAMPLE  VALUE 

DESCRIPTION 

c 

3.0  x10s  ms'1 

speed  of  light 

Re 

x  6.37  x  1 06  m 

effective  earth  radius 

Mc 

0.1  g  m'3 

cloud  liquid  water  content 

Kf 

0.92 

complex  refractive  index  function  for  water 

a 

200  mm' m'3  (mm  h'1)'6 

Z-R  power  law  coefficient* 

b 

1.6 

Z-R  power  law  exponent* 

h0 

4000  m 

height  of  0  °C  level 

hm 

20000  m 

maximum  height  of  precipitation 

hb 

25  m 

cloud  base  height 

h 

600  m 

thickness  of  the  melting  layer 

Z(h ) 

see  Figure  3 

reflectivity  factor  profile 

TSea 

28  °C 

sea  surface  temperature 

S 

34ppt 

sea  surface  salinity  in  parts  per  thousand  (ppt) 

Hs 

1 

sea  state 

'Z=aRb 
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Reflectivity  Factor  Profile 


Figure  3.  Example  height  profile  of  reflectivity  factor,  Z,  where  dBZ  =  10.0  log(Z  in  mm6  m  3) 


4.0  km  is  where  the  snow  melts  into  rain.  The  rain  only  region  is  below  3.4  km  altitude. 

The  main  effect  of  the  precipitation  field  on  the  target  signal  is  attenuation  of  the  received 
power.  The  amount  of  attenuation  is  a  function  of  the  transmitted  frequency  and  polarization, 
the  cloud  liquid  water  content,  the  precipitation  rates,  and  the  specific  path  travelled  by  the 
radar  signal  through  the  atmosphere.  The  end  of  the  signal  path  is  the  location  of  the  target  as 
specified  by  the  range  vector  rt .  Thus,  an  attenuation  value,  La  (rt ) ,  is  calculated  for  each 
missile  position  or  equivalently  for  each  transmitted  radar  pulse.  The  empirical  models  and 
mathematical  approach  used  for  these  calculations  are  described  in  detail  in  Thomson  [1]. 
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4  Radar  Model 


The  radar  model  of  this  simulation  is  very  similar  to  that  which  was  described  for  the  weather 
clutter  simulation  in  Thomson  [1].  Most  of  the  differences  are  a  result  of  upgrades  that  were 
added  to  allow  the  radar  model  to  be  more  generally  applicable.  Upgrades  related  to  the  range 
weighting  function  and  the  antenna  gain  pattern  have  also  been  implemented  within  the 
current  version  of  the  weather  clutter  simulation.  Thus,  the  descriptions  given  in  Sections  4.3 
and  4.4  also  apply  to  the  current  version  of  the  weather  clutter  simulation. 

The  main  remaining  difference  between  the  two  radar  models  is  the  complexity  of  the 
transmitted  waveform.  The  model  presented  here  considers  a  multiple  burst  waveform  in 
which  the  frequency,  pulse  repetition  time,  and  the  number  of  pulses  per  coherent  burst 
change  from  burst  to  burst.  The  bursts  are  grouped  into  dwells  so  that  an  M-of-N  detection 
scheme  can  be  applied  later.  The  waveform  parameters  describing  the  bursts  can  also  vary 
from  dwell  to  dwell.  In  addition,  all  of  the  signal  samples  for  all  dwells  are  calculated 
together  for  the  entire  target  approach.  In  the  case  of  the  weather  clutter  simulation  only  a 
single  burst  is  simulated.  The  more  complex  waveform  used  in  this  target  simulation  can  still, 
however,  be  modelled  by  the  weather  clutter  simulation,  but  this  requires  the  weather  clutter 
simulation,  as  it  has  been  described  in  Thomson  [1],  to  be  executed  multiple  times;  once  for 
each  burst.  Thus,  the  outputs  of  the  two  simulations  are  compatible  and  can  be  combined  as 
described  in  Thomson  [3]. 

4.1  System  parameters 

The  parameters  used  to  describe  a  radar  system  are  shown  in  Table  2.  These  example 
waveform  parameters  only  specify  a  single  burst.  Section  4.2  discusses  the  full  flexibility  of 
the  transmitted  waveform.  The  peak  gain  on  boresight,  as  well  as  the  azimuth  and  elevation 
beamwidths,  are  calculated  as  a  function  of  frequency  as  discussed  in  Section  4.4. 

Pulse  compression  is  modelled  using  a  simplistic  approach  (see  Barton  [7],  p  367)  similar  to 
that  used  in  the  weather  clutter  simulation  (Thomson  [1]).  The  compressed  pulse  is  assumed 
to  have  the  same  properties  as  the  uncompressed  pulse,  except  for  its  length  and  transmitted 
power.  The  compressed  pulse  length  is  given  by  rc  =  r/Pc ,  where  Pc  is  the  pulse 
compression  ratio  and  ris  the  uncompressed  pulse  length.  The  compressed  pulse  width  is 
then  used  in  all  calculations,  such  as  the  range  weighting  function.  This  approach  ignores  the 
effects  of  range  sidelobes.  The  increase  in  transmitted  power  associated  with  pulse 
compression  is  modelled  by  including  a  factor  of  Pc  in  the  equation  used  to  calculate  received 
power  (see  Equation  (36)).  This  approach  assumes  that  the  missile  is  a  coherent  target  for  a 
time  period  greater  than  r. 

In  the  case  of  the  weather  clutter  simulation  as  described  in  Thomson  [1],  pulse  compression 
was  implicitly  modelled  by  inputting  appropriately  scaled  parameters  into  the  simulation,  i.e., 

r/Pc  ,  was  entered  for  the  pulse  length,  and  Pc2Pt  was  entered  for  the  transmitted  power.  In 
this  simulation,  however,  an  input  parameter  for  pulse  compression  is  required  and  the 


8 


DRDC  Ottawa  TR  2002-145 


Table  2.  Parameters  describing  the  radar. 


SYMBOL 

EXAMPLE  VALUE 

DESCRIPTION 

ha 

25  m 

antenna  height  above  sea  level 

Vc 

15° 

antenna  cant  back  angle 

0b 

100° 

azimuth  of  antenna  boresight 

Pt 

17.12  kW 

transmitted  power 

0 

90°  (vertical) 

transmitted  polarization 

T 

50  ns 

uncompressed  pulse  width 

Pc 

250 

pulse  compression  ratio 

Tc 

200  ns 

compressed  pulse  width 

60  dB 

receiver  system  gain 

b6 

4  MHz 

6  dB  receiver  bandwidth 

G't 

37.45  dB 

antenna  gain  on  boresight  at  fg 

f. 

10.0  GHz 

frequency  corresponding  to  boresight  gain  specification 

Gf 

1.0459076 

fit  parameter  for  antenna  gain  pattern 

06b  {ft ) 

2.6471614 

6  dB  azimuth  beamwidth  on  boresight 

06b  (ft) 

2.2864836 

6  dB  elevation  beamwidth  on  boresight 

A  0s 

90° 

azimuth  extent  of  scan  area 

ec 

0° 

elevation  of  beam  axis  when  target  is  illuminated 

(frc 

130° 

azimuth  of  beam  axis  when  target  is  illuminated 

r min 

7.545  km 

minimum  range  of  horizon  search 

r max 

29.745  km 

maximum  range  of  horizon  search 

A  r 

30  m 

distance  between  the  centres  of  adjacent  range  bins 

Ts 

240.7  ns 

pulse  repetition  time 

ft 

10.93  GHz 

transmitted  frequency 

Nc 

26 

number  of  pulses  per  coherent  burst 

necessary  calculations  are  performed  within  the  simulation.  Note  that  in  Thomson  [1]  (see 
Section  6.1)  Pt  was  scaled  by  Pc2  whereas  in  this  document  Pt  is  scaled  by  Pc  to  model 
pulse  compression.  This  scaling  ( Pc 2 )  could  also  be  considered  to  give  correct  results  if  the 
extra  factor  of  Pc  is  thought  of  as  system  gain,  since  in  Thomson  [1]  the  noise  figure  was  also 
scaled  by  a  factor  of  Pc .  However,  this  extra  scaling  is  unnecessary  and  is  no  longer  used. 
Pulse  compression  in  the  upgraded  weather  clutter  simulation  would  now  be  implemented  by 
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using  the  compressed  pulse  length  as  the  input  pulse  length  and  PcPt  as  the  input  transmitted 
power.  The  noise  figure  would  not  be  scaled. 

4.2  Transmitted  waveform 

The  transmitted  waveform  is  defined  by  the  parameters  shown  in  Table  2.  However,  as  stated 
earlier  T„f,  and  Nc  can  be  varied  from  burst  to  burst  and  from  dwell  to  dwell.  Three  levels  of 
complexity  are  possible.  First  the  waveform  parameters  can  remain  constant  so  that  only  a 
single  burst  is  defined  and  is  transmitted  repeatedly  until  the  target  reaches  the  minimum 
range,  rmin.  Second,  a  set  of  Nb  values  for  each  of  Ts,f,  and  Nc  can  be  defined.  This  results  in 
a  transmitted  dwell  consisting  of  Nb  bursts.  This  dwell  waveform  is  then  transmitted 
repeatedly.  Finally,  a  set  of  Nb  x  Nj  values  can  be  defined  for  each  of  Ts,f,  and  Nc.  This 
results  in  Nd  dwells  being  transmitted  where  each  dwell  contains  Nb  bursts  and  all  bursts  can 
be  unique.  If,  in  this  case,  the  required  transmission  time  for  all  defined  dwells  is  shorter  than 
the  time  required  for  the  target  to  reach  the  minimum  range,  then  the  dwells  will  be  repeated 
in  order  of  their  definition. 


4.3  Range  weighting  function 

An  individual  I  or  Q  sample  corresponding  to  a  given  range  bin  represents  the  integration  of 
the  received  signal  voltage  over  a  finite  period  of  time,  which  corresponds  to  a  finite  range 
increment.  Therefore,  the  backscattered  signals  from  all  scatterers  within  this  range  increment 
will  contribute  to  the  signal  sample.  The  range  weighting  function,  W2,  describes  the  relative 
contribution  to  the  signal  sample  of  each  scatterer’s  backscattered  power  based  on  the  position 
of  the  scatterer  with  respect  to  the  range,  rc,  corresponding  to  the  sample  time. 

The  only  scatterers  considered  in  this  simulation  are  the  different  scattering  centres  of  the 
single  target  approaching  the  radar.  These  scatterers  are  modelled  as  a  single  point  scatterer 
with  the  effect  of  the  multiple  scattering  centres  of  the  target  being  generated  statistically  (see 
Section  5).  Thus,  the  entire  target  is  assumed  to  be  completely  within  the  resolution  volume 
corresponding  to  the  signal  sample.  The  weight,  W2(rh  rc),  for  the  point  target  at  range,  rt,  is 
applied  to  the  received  power  calculated  via  the  radar  equation  (see  Equation  (36)). 


The  specific  form  of  the  range  weighting  function  used  in  this  simulation  is  taken  from 
Doviak  and  Zmid  [6],  i.e., 


W2(/-,rc)=^jerf 


c-Jh\2 


r„  -r  +  - 


CTr 


-erf 


V 


\c 


Vln2 


r-r-- 


ct • 


A 


(7) 


where  c  is  the  speed  of  light,  and  “erf”  is  the  error  function  given  by 


o 


(8) 
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Equation  (7)  corresponds  to  the  case  of  a  rectangular  transmitted  pulse  and  a  practically 
matched  receiver  having  a  Gaussian  frequency  response. 

This  form  of  the  range  weighting  function  is  also  used  in  the  current  version  of  the  weather 
clutter  simulation.  The  single  difference  between  Equation  (7)  and  its  counterpart  that  is 
described  in  Thomson  [1]  is  the  decoupling  of  B6  and  Tc. 

Figure  4  shows  a  plot  of  Equation  (7)  for  the  case  defined  by  the  parameters  in  Table  2.  The 
6  dB  range  width  of  the  resolution  volume,  is  determined  by  solving  the  transcendental 
equation 


W‘ 


rcr-rc 


(9) 


For  B6  =  4.0  MHz,  rc  =  200  ns,  and  any  value  of  rc,  Rf,  =  39.806892  m,  which,  as  shown  in 
Figure  4,  is  greater  than  cTc/ 2  =  30  m  .  The  target  simulation  only  considers  the  portion  of 
the  resolution  volume  within  the  6  dB  range  points,  rc±R6/ 2 ,  as  is  the  case  for  the  weather 
clutter  simulation.  For  the  case  defined  in  Table  2,  R^  is  also  greater  than  the  distance 
between  the  centres  of  adjacent  range  bins,  A r.  Thus,  in  this  case,  there  is  some  overlap 
between  adjacent  range  bins. 


4.4  Antenna  gain  pattern 

The  radar  antenna  is  modelled  as  an  elliptical  phased  array.  The  corresponding  radar  beam 
has  frequency  dependent  beamwidths  that  increase  and  a  frequency  dependent  peak  gain  that 


Figure  4.  The  range  weighting  function. 
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decreases  as  the  beam  is  scanned  away  from  the  boresight. 

The  geometry  of  the  phased  array  antenna  is  shown  in  Figure  5.  Here  the  antenna  face  is 
shown  as  an  elliptical  aperture  canted  back  by  an  angle  \ffc  with  respect  to  the  z-axis.  The  unit 
vector  rb  defines  the  direction  of  the  antenna  boresight  in  a  coordinate  system  that  has  its 

origin  at  the  centre  of  the  antenna.  The  angle  specifies  the  azimuth  of  the  antenna 
boresight. 

The  target  simulation  only  considers  the  portion  of  the  resolution  volume  that  receives  a 
weight,  specified  by  the  radar  transmit/receive  characteristics,  that  is  within  6  dB  of  the 
weight  applied  to  the  centre  of  the  resolution  volume  at  spherical  polar  coordinates  ( rc ,  6C,  <p,). 
Therefore,  the  antenna  gain  pattern  must  be  specified  over  the  angular  range  defined  by  the 
3  dB  elevation  and  azimuth  beamwidths,  $}  and  A  two  dimensional  Gaussian  function  is 
used  to  describe  the  far-field  1-way  antenna  gain  pattern  in  this  region,  i.e., 


z 


Figure  5.  Antenna  geometry. 


12 


DRDC  Ottawa  TR  2002-145 


(10) 


G{6, <p,0c,<pc)  =  Gf  Gs exp 


-81n2 


Q-Oc\2  l^-^fsin2^ 

%  fe 


where  #and  0are  the  polar  and  azimuth  angles,  respectively,  that  specify  the  location  of  the 
point  in  the  spherical  polar  coordinate  system  for  which  the  gain  is  to  be  calculated.  The  polar 
and  azimuth  angles  specifying  the  radar  beam  axis  are  given  by  9C  and  <pc,  respectively.  The 

6  dB  beamwidths,  06  =  sj2&3  and  (f>6  =  yf203 ,  specify  an  angular  distance  in  a  coordinate 

system  that  is  fixed  with  respect  to  the  radar  beam  axis.  The  quantities  1 9  -  0C\  and  |  </>  -  </>c\ , 

however,  specify  angular  distances  in  the  spherical  polar  coordinate  system.  Thus,  the  sin2  9 
term  in  Equation  (10)  is  required  to  convert  azimuthal  distances  in  the  spherical  polar 
coordinate  system  to  angular  distances  relative  to  the  antenna  gain  pattern.  The  dimensionless 
parameter  G/is  included  to  improve  the  fit  of  this  Gaussian  model  to  real  antenna  patterns. 

The  peak  gain  on  the  radar  beam  axis  is  given  by  Gs.  The  decrease  in  peak  gain  caused  by 
scanning  the  beam  off  boresight  is  given  by  Billetter  [8]  as 

Gs=Gb(cosycf5  (11) 

where  Gb  is  the  gain  along  the  antenna  boresight  and  yc  =  cos-1  (re  •  rb )  is  the  scan  angle.  The 
unit  vector  rc  specifies  the  direction  of  the  centre  of  the  resolution  volume,  hence  yc  is  given 
by 


yc  =  cos  1  (sin 9C  sin <pc  cos  y/c  sin^,  +  sin 9C  cos  </>c  cos  y/c  cos  <ph  +  cos#c  sin y/c)  (12) 


Thus,  the  loss  in  peak  gain  is  a  function  of  the  magnitude  of  the  scan  angle  and  not  a  function 
of  the  direction  of  this  scan  angle,  even  though  the  antenna  array  is  elliptical.  This  approach 
implicitly  assumes  that  the  radiation  patterns  of  the  individual  elements  that  make  up  the 
antenna  array  are  circularly  symmetric. 

The  broadening  of  the  beam  due  to  scanning,  however,  is  a  function  of  effective  aperture 
reduction  (see  Billetter  [8]),  which  is  direction  dependent.  Thus,  the  increases  in  61,  and 
caused  by  scanning  off  boresight  are  calculated  separately  as 


*6  = 


6b 


n 


\ 


cos  i  ~-vc-e, 

2  J 


“6b 


cos  {<t>b  -<t>c) 


(13) 


where  the  6  dB  beamwidths  on  boresight,  06 b  and  ((kb,  are  required  input  parameters  that  are 
specified  as  single  values  or  as  sets  of  frequency  dependent  values.  In  the  case  where  a 
frequency  dependent  set  is  used,  (e.g.,  { 06b}  =  {3.2°,  3.0°,  2.8°,  2.6°,  2.4°,  2.2°,  2.0°}  for 
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{ft}  =  {8.5  GHz,  9.0  GHz,  9.5  GHz,  10.0  GHz,  10.5  GHz,  11.0  GHz,  11.5  GHz})  the 
beamwidth  values  for  a  specific  transmitted  frequency  are  calculated  by  linear  interpolation 
within  the  sets. 

The  gain  along  the  antenna  boresight,  Gb,  is  also  dependent  upon  the  transmitted  frequency,/, 
(see  Doviak  and  Zmid  [6]),  i.e., 


G 


4  nAef? 


(14) 


where  Ae  is  the  effective  area  of  the  antenna.  The  input  parameters  G'b  and/  (see  Table  2), 
which  specify  the  boresight  gain  for  a  particular  frequency,  can  be  used  to  calculate  Ae  via 
Equation  (14).  Thus,  the  boresight  gain  for  arbitrary  transmitted  frequencies  is  calculated  as 

c»(/,)=c;4  d5) 

Jg 

As  was  the  case  for  the  range  weighting  function,  the  current  version  of  the  weather  clutter 
simulation  has  been  upgraded  from  that  described  in  Thomson  [1]  to  use  this  model  for  the 
antenna  gain  pattern. 


4.5  Scan  pattern 

The  modelled  radar  electronically  scans  in  azimuth  only.  The  azimuth  extent  of  the  scanned 
area  is  given  by  A<ps .  Figure  6  depicts  how  this  azimuth  interval  is  observed  by  a  set  of 
overlapping  radar  beams  that  adhere  to  the  phased  array  model  outlined  in  Section  4.4. 

The  number  of  required  beams,  A^ram,  is  determined  by  dividing  A<ps  with  the  smallest  3  dB 
beamwidth  on  boresight  that  is  associated  with  the  defined  waveform  set,  i.e., 


N 


beam 


(16) 


where  int  takes  the  integer  of  its  argument,  min  takes  the  minimum  of  its  argument,  and 
{&  J  =  {U/V2  represents  the  set  of  3  dB  beamwidths. 


The  azimuth  positions  of  the  beams  are  set  by  spacing  their  beam  axes  apart  by  min{  fob). 
The  beams  are  transmitted  in  order  of  increasing  beam  axis  azimuth.  Calculation  of  the 
received  signal  samples,  however,  is  only  performed  for  the  particular  beam  in  which  the 
target  is  approaching  (e.g.,  the  yellow  shaded  beam  in  Figure  6).  The  other  beams  are 
considered  only  for  determining  the  time  delay  between  measurements  along  the  beam 
containing  the  target  (see  Section  5.1). 
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Figure  6.  Azimuthal  scan  pattern.  The  depicted  3  dB  beamwidth  on  boresight  is  10°.  The  azimuth 
extent  of  the  scan  area  is  90°.  The  dashed  lines  show  the  boundaries  of  the  scan  area  and  the  antenna 
azimuth  boresight.  Grey  shaded  regions  depict  the  individual  beams.  The  beam  for  which  radar  signals 

are  calculated  is  shaded  yellow. 
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5  Missile  Model 


5.1  Trajectory 

The  modelled  missile  flies  at  a  constant  height  above  sea  level  and  at  a  constant  velocity 
directly  toward  the  radar.  The  input  parameters  required  to  define  the  missile  properties  are 
given  in  Table  3.  The  trajectory  of  the  missile  is  defined  in  terms  of  a  set  of  points, 
rt  =  (rf,  0t,  fa),  that  define  the  location  of  the  missile  each  time  it  is  illuminated  by  a  radar 
pulse.  To  calculate  these  points  the  ranges,  a,,  along  the  missile  trajectory  are  first 
determined. 

Figure  7  depicts  the  geometry  of  the  missile  trajectory.  The  radar  antenna  is  located  at  the 
origin  of  the  y  and  z  axes.  The  starting  range  of  the  missile,  rs,  is  a  required  input  parameter 
of  the  target  simulation.  The  distance  along  the  missile  trajectory  from  the  starting  point  to 
the  point  where  the  missile  trajectory  intercepts  the  z  axis  is  given  by 

as={Re  +  h,)/3s  (17) 

where  Re  is  the  effective  earth  radius  and  h,  is  the  height  of  the  missile  above  sea  level.  The 
angle  fis  can  be  related  to  rs  using  the  law  of  cosines,  i.e., 

t  =  (K  +h,f+  (*,  +Kf-  2{Re  +  h,  \Re  +  ha  )cos  J3S  (18) 

where  ha  is  the  height  of  the  radar  antenna  above  sea  level.  Thus,  as  can  be  written  in  terms  of 
known  quantities  as 

as={Re  +  h,)cos,~l 

The  scenario  is  considered  to  end  when  the  missile  reaches  the  range,  rec  =cr/2 ,  where 
eclipsing  begins.  Here  r  is  the  uncompressed  pulse  width.  The  range  along  the  flight  path 


Table  3.  Parameters  describing  the  target. 


SYMBOL 

EXAMPLE  VALUE 

DESCRIPTION 

N/A 

3 

Swerling  target  type 

<*b 

0.01  m2 

backscatter  cross-section 

v, 

281  m  s'1 

target  speed 

h, 

20  m 

height  of  target  above  sea  level 

0s 

131.17324° 

azimuth  of  target  trajectory 

rs 

30901.541  m 

target  start  range 

{Re+ht)2+{Re+ha)2-r? 

2{Re+h,){Re+K) 


(19) 
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Figure  7.  Geometry  of  the  missile  trajectory.  The  radar  antenna  is  located  at  the  origin  of  the  y  and  z 
axes.  The  thick  blue  curve  represents  the  surface  of  the  sea. 


from  the  z  axis  to  rec  is  calculated  similarly  to  as,  i.e., 


a  =(R  +h)  cos-1  ^Re  +fy).  +  K) _ r_ec_ 

aec  {Re  +  ht)  y  2{Re+ht){Re+ha)  J 


The  first  radar  pulse  of  the  scenario  is  transmitted  along  azimuth  (f)s  when  the  missile  is  at 
position  rs  or  at  range  as  along  the  flight  path.  Therefore,  the  number  of  complete  scans  over 
the  azimuth  interval  A</>s  that  the  radar  will  perform  during  the  scenario,  or  equivalently  the 
number  of  transmitted  dwells,  Nj,  is  given  by 


where  v,  is  the  missile  speed,  p  is  an  index  that  is  used  to  sum  over  the  individual  dwells  of 
the  scenario,  q  is  an  index  that  is  used  to  sum  over  the  individual  bursts  of  a  dwell,  Tspq  is  the 
pulse  repetition  time  for  the  qh  burst  of  the  pth  dwell,  and  Ncpq  is  the  number  of  coherent 
pulses  in  the  qth  burst  of  the  plh  dwell.  Equation  (21)  indicates  that  each  dwell  waveform  is 
transmitted  along  all  azimuths  until  the  scan  is  complete,  at  which  time  the  dwell  waveform  is 
changed  to  the  next  defined  waveform. 
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The  transmission  time  of  the  first  pulse  of  the  scenario  is  set  to  t0oo  =  0.0  s.  The  processing 
power  of  the  radar  system  is  assumed  to  be  great  enough  such  that  there  are  no  delays 
between  transmitted  bursts  or  dwells,  i.e.,  the  transmission  times  are  completely  determined 
by  the  defined  nc  and  Ts  values.  Thus,  the  transmission  times  of  all  subsequent  pulses  are 
calculated  as 


lkmn 


■■N 


k  Ny- 1 

beamY,LV- 

p= 0  q=0 


N 

spq  cpq 


skin  ^  ckm  'skm 


9=0 


9=0 


n=0,\,...,Nckm  -1 
m  =  0, 1  ,...,Nb  -1 
k=0,l,...,Nd-l 

Here  thm  is  the  transmission  time  for  the  pulse  of  the  m*  burst  of  the  dwell.  The  missile 
positions  along  the  flight  path  at  the  time  each  pulse  is  transmitted,  a',(thnn ) ,  are  then 
calculated  as 


a,t(tkmn)=as-\ltkmn  (23) 

The  transmitted  pulses  will  intercept  the  missile  at  ranges  )  such  that 

( t/ann  )  “  (  ^knm  )  —  ^  f  ^kinn  (24) 

where  Atbnn  are  the  times  required  for  the  transmitted  pulses  to  travel  from  the  antenna  to  the 
missile  at  ranges  ).  Using  the  law  of  cosines,  Equation  (24),  and  /?,  -at/{Re  +  ht), 

rt{tkmn)  can  be  written  as 


r, ( ) =  («e  +  h,  y 2  +  (Re  +  ha  y 2  -  2(Re  +  h,  ){Re  +  ha )  cos 


a’Atbm)-v1 At, 


kmn 


R„  +  h, 


(25) 


The  ranges  to  the  missile  positions,  as  defined  by  Equation  (25),  are  also  equal  to  cAthm  , 
where  c  is  the  speed  of  light.  Thus,  the  remaining  unknowns  in  Equation  (25),  Athnn  ,  are 
determined  by  numerically  solving  the  transcendental  equations  that  are  formed  by  equating 
the  right  hand  side  of  Equation  (25)  to  cAt^.  Substituting  the  resultant  Athnn  values  back 
into  Equation  (25)  will  allow  the  ranges  required  for  calculating  the  received  power  (see 
Section  5.2)  to  be  determined. 

Once  the  r((  )  have  been  calculated  the  law  of  cosines  is  again  used  to  calculate  the  polar 

angles  of  the  missile  positions  at  the  times  that  the  missile  is  illuminated,  i.e., 
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5.2  I  and  Q  signal  samples 

The  radar  pulses  of  the  transmitted  waveform  have  been  grouped  into  bursts  and  dwells.  The 
missile  positions  corresponding  to  these  pulses  are  also  grouped  together  according  to  the 
bursts  and  dwells  of  the  transmitted  waveform.  Depending  upon  its  speed  and  initial  position, 
the  missile  may  move  across  more  than  one  range  bin  during  the  time  required  to  transmit  a 
burst.  In  addition,  the  measurement  ray  may  be  defined  such  that  adjacent  range  bins  overlap. 
Thus,  for  each  burst,  time  series  of  I  and  Q  signal  samples  must  be  calculated  for  all  range 
bins  that  contain  at  least  one  missile  position  associated  with  the  burst. 

Figure  8  depicts  an  example  scenario  in  which  the  missile  positions  of  the  burst  affect  more 
than  one  range  bin.  In  this  example  the  adjacent  range  bins  overlap  and  the  missile  is  moving 
fast  enough  to  move  across  two  range  bins.  The  mth  burst  of  the  Jr  dwell  has  only  four  radar 
pulses  (Nchn  =  4)  and  consequently  only  four  corresponding  missile  positions.  These  four 
positions  span  across  range  bins  A  and  B  (shaded  yellow),  which  are  defined  by  their  6  dB 
resolution  volumes.  Thus,  two  sets  of  I  and  Q  time  series  will  be  calculated  for  this  burst. 
However,  since  the  calculations  only  consider  contributions  to  the  signal  samples  from  targets 
within  the  6  dB  resolution  volume,  I4  =  Q4  =0  for  the  time  series  corresponding  to  range 
bin  A,  and  h  =  =  I2  =  Q2  =  0  for  the  time  series  corresponding  to  range  bin  B.  Therefore, 

the  average  power  received  from  the  burst  is  likely  to  be  greater  for  the  time  series  of  range 
bin  A  than  that  for  range  bin  B  even  though  range  bin  A  is  further  in  range. 

The  components  of  an  individual  complex  signal  sample  corresponding  to  a  given  missile 
position  and  a  given  range  bin,  are  calculated  as 


I  =  -Jf\  cos  rp 
Q  =  J1\  sirup 


(27) 


Figure  8.  Example  of  burst  missile  positions  relative  to  the  radar  resolution  volumes. 
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where  Pr  is  the  power  received  from  the  missile,  and  pis  the  phase  difference  between  the 
radar  reference  signal  and  the  received  voltage  sample.  For  simplicity,  the  subscripts 
denoting  the  dwell,  burst,  and  pulse  have  been  omitted.  Equation  (27)  assumes  that  the 
resistance  of  the  radar  receiver  is  1  Q. 


The  phase  of  the  reference  signal  at  time  Amo  =  0-0  s  is  assumed  to  be  equal  to  zero.  Thus,  the 
value  of  (p  that  represents  the  relative  phase  of  a  particular  I  and  Q  pair  is  calculated  as 


—  Aw.  r 


(28) 


where  r,  is  the  corresponding  missile  range,  and  A  is  the  transmitted  wavelength  of  the  burst 
for  which  the  I  and  Q  time  series  are  being  calculated. 

The  calculation  of  the  received  power,  P,,  corresponding  to  a  particular  I  and  Q  pair,  is  based 
upon  the  integral  form  of  the  radar  equation  that  was  used  in  Thomson  [1]  (see  also  Doviak 
and  Zmid  [6]).  The  power  received  from  a  collection  of  scatterers  within  the  range  bin 
centred  at  rc  =  (rc ,  0C ,  <f>c )  is  determined  by  integrating  the  weighted  reflectivity  field  over  the 
6  dB  resolution  volume,  v,  i.e., 


Pr{*c)=\v{T)Y{T,?c)  d»  (29) 

V 

where  T]  describes  the  reflectivity  field,  and  Y  is  a  composite  weighting  function  dependent 
upon  the  transmit/receive  characteristics  of  the  radar  system.  In  the  target  simulation, 
however,  there  is  only  one  point  scatterer  located  at  position  rf .  Thus,  the  reflectivity  field  is 
given  by 

ifc)’oAr -nha" 00) 

rt  sin  0, 

where  S(r  -  rt)  is  the  Dirac  delta  function  defined  by  (see  Jordan  and  Balmain  [9]) 


j  S(r  -  rt)dv  =j^ 


if  rt  is  in  v 
if  rt  is  not  in  v 


(31) 


The  backscatter  cross  section  of  the  target,  Oj,  is  an  input  parameter  required  by  the 
simulation.  However,  its  value  is  set  to  zero  for  all  target  positions  that  are  below  the  radar 
horizon,  i.e.,  for  target  positions  such  that 
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6t  >  k/2 


(Re  +  ha )  sin(;r  -  et)<Re 


rt>£s 

where  £s  is  the  distance  from  the  radar  antenna  to  the  sea  surface  along  the  path  through  the 
earth  to  the  missile  position  (see  Figure  9).  By  applying  the  law  of  cosines  to  the  left  triangle 
shown  in  Figure  9,  ls  can  be  written  as 

ls  =[Rl+{Re  +  hay -2Re{Re+ha)sm(et-¥g)Y2  (33) 

The  grazing  angle  y/g is  calculated  by  applying  the  law  of  sines  to  the  left  triangle  of  Figure  9, 
i.e., 


^=--sm 


7t  .  _i  (  R  e  +  ha  |  .  /  \ 

--sm  sm{7t-9t) 


Substituting  Equation  (30)  into  Equation  (29)  gives 


pr  0 ) : =  J 1 a b  5  (r  -  rt ) : y  ( i r ,  rc )  ■ d  v  =  a bY  {ft,  rc ) 


z 


Figure  9.  Geometry  describing  target  positions  below  the  radar  horizon.  The  thick  blue  curve 

represents  the  sea  surface. 
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when  rt  is  within  v.  Expansion  of  the  composite  weighting  function  (Thomson  [1], 

Doviak  and  Zmic  [6])  gives  a  form  of  the  radar  equation  for  a  point  target  that  allows  the 
position  of  the  target  within  the  resolution  volume,  and  its  corresponding  weighting  by  the 
antenna  gain  pattern  and  the  radar  receiver,  to  be  considered,  i.e., 

PCP,  #  Gsy,  G1  (&,,</>,, 8c,<pc)  w2(r„rc)ab 

Prvc)  = -  (TxTn  i".  x- -  (36> 

M  LaVn0n<P,)r, 

where  P,  is  the  transmitted  power,  and  Gsys  is  the  receiver  system  gain.  Note  that  a  term,  Lr, 
was  incorrectly  included  in  the  composite  weighting  function  of  Thomson  [1]  to  describe  a 
power  loss  due  to  finite  receiver  bandwidth.  This  term  is  not  necessary  when  the  range 
weighting  function  is  included  in  the  composite  weighting  function.  Hence  Lr  does  not  appear 
in  Equation  (36)  and  is  no  longer  a  parameter  of  the  weather  clutter  simulation  (Thomson  [1]). 
The  factor  Pc  has  been  included  to  model  the  increase  in  power  backscattered  from  a  coherent 
target  when  pulse  compression  is  used.  The  2-way  attenuation  of  the  radar  signal  over  the 
path  from  the  antenna  to  the  target  and  back,  La  (r( ,  0t ,  0, ) ,  is  calculated  using  the  empirical 
formulas  described  in  Thomson  [1]. 

5.3  Multipath  interference 

Refraction  of  the  radar  signals  in  the  atmosphere  has  been  accounted  for  by  using  an  effective 
earth  radius  in  all  calculations.  This  approach  allows  the  signal  ray  paths  to  be  considered  as 
straight  lines,  which  simplifies  the  calculations  required  to  determine  the  interference  effects 
caused  by  a  reflective  sea  surface.  In  this  case,  the  radar  signal  will  travel  along  four  different 
paths  to  the  target  and  back  to  the  antenna.  These  paths  are  shown  in  Figure  10. 

The  calculations  described  in  Section  5.2  apply  to  the  signals  received  along  path  1.  The  I 
and  Q  signals  samples  corresponding  to  paths  2,  3,  and  4  are  calculated  in  a  similar  way, 
except  for  adjustments  related  to  the  reflection  coefficient,  the  path  length,  the  phase  shift 
upon  reflection,  and  the  divergence  of  the  radar  beam  caused  by  reflection  from  a  spherical 
surface. 


Path  3  involves  a  reflection  of  the  radar  signal  from  the  sea  surface  both  on  the  way  to  the 
target  and  on  the  way  back.  For  both  reflections  the  signal  is  reflected  from  a  medium  (sea) 
that  is  more  optically  dense  than  the  medium  in  which  it  is  travelling  (air).  Therefore,  each 


Figure  10.  The  signal  paths  that  cause  interference.  The  thick  blue  curve  represents  the  sea  surface. 
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reflection  will  result  in  a  phase  shift  of  n  radians.  Consequently,  the  phase  of  the  received 
signal  can  be  calculated  as  in  Equation  (28),  i.e., 


<p3  = 


(37) 


where  rt3  =  it3i  +  lt3r  is  the  range  along  path  3.  Figure  1 1  defines  the  quantities  £m  and  it3r 
as  the  range  from  the  antenna  to  the  reflection  point  and  the  range  from  the  reflection  point  to 
the  target,  respectively. 

The  power  received  from  the  target  along  path  3  is  calculated  using  the  radar  equation,  i.e., 


PrArc)=: 


p2PcPt  A2  Gsys  G2{9t3.^dc,(Pc)  w2(rl3,rc)crb 


(4/r)34(rf3,0,3,$)rr43 


(38) 


However,  in  this  case  the  received  power  will  be  reduced  based  upon  the  reflectance  of  the 
sea  surface,  p.  In  addition,  the  composite  weighting  function  must  be  evaluated  at  the  polar 
angle  that  specifies  the  direction  of  the  reflection  point,  dt3 ,  and  at  the  range  rf3 . 


The  calculation  of  p,  £m ,  lt3r ,  and  dl3  in  the  target  simulation  is  done  in  a  very  similar  way 
to  that  described  for  the  weather  clutter  simulation  in  Thomson  [1].  The  main  difference 
between  the  two  calculations  is  the  divergence  factor.  These  calculations  are  reproduced  in 
Annex  A  of  this  report  because  of  this  difference  and  also  because  of  differences  in  notation. 


The  I  and  Q  components  corresponding  to  signals  received  along  path  3  are  calculated  as 


I3n=V^TCOS^3n 

Q3»  =4P^Sin(P3n 


(39) 


In  Equation  (39),  the  subscript  “3”  refers  to  path  3  and  the  subscript  refers  to  the 
individual  signal  samples  of  a  received  burst  as  in  Section  2. 

Signal  path  2  involves  only  one  reflection  from  the  sea  surface.  Thus,  a  phase  shift  of  7t 
radians  must  be  added  to  the  phase  calculation,  i.e., 


<p2=x 


A 


(40) 


where  the  apparent  range  is  calculated  as  rt2  =  (rt  +  rl3)/2.  The  single  reflection  is  also 
represented  in  the  power  calculation  as  a  factor  of  p,  i.e., 
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Figure  11.  Path  3  geometry. 


„  ,^_pr.r,?omGl  iv’fe.'iK 

rr2  Vc ) - 77~v5 - 7 - 

(4/v)  Lal  rt2 

In  this  case  the  apparent  gain  corresponding  to  path  2,  G2,  is  calculated  as 

G2=MO'>tn0c'tc)  G{0t3,tt,ffc,tc) 

and  the  apparent  2-way  attenuation  corresponding  to  path  2,  La2 ,  is  calculated  as 


(41) 


(42) 
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(43) 


Ai2  — 


Lg(rtA’</>t)+  Lg(rtlA3’<P,) 

2 


The  I  and  Q  samples  corresponding  to  signal  path  4  are  exactly  the  same  as  those  for  signal 
path  2.  Hence,  the  path  2  and  4  signal  samples  are  given  by 


I4 n  ~  I In  -  *JPr2n  C0S  ^2n 

Q4n=Q2n=4P^Sia<P2n 


(44) 


Once  the  signal  samples  have  been  calculated  for  the  four  different  paths,  the  overall  received 
signal  samples  for  a  given  burst  and  a  given  range  bin  are  calculated  by  coherent  summation, 
i.e., 


~  Iln  +  ^2n  +  ^3n  +  ^4 n 
Qn  =  Ql  n  +  Q2  n  +  Q3n  +  Q4n 


(45) 


As  shown  in  Figure  8,  the  received  pulses  corresponding  to  a  particular  burst  can  contribute 
backscattered  power  to  the  signal  samples  corresponding  to  more  than  one  range  bin.  In 
general,  it  is  possible  that  the  set  of  range  bins  associated  with  a  given  burst  may  not  be  the 
same  for  all  of  the  different  signal  paths  because  of  their  differing  path  lengths.  The  target 
simulation  properly  considers  these  range  differences  when  summing  the  multipath  signal 
samples  for  each  range  bin. 
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6  Statistical  Model 


As  shown  in  Figure  2  the  statistical  model  creates  random  numbers  that  will  be  filtered  by  the 
output  of  the  missile  model.  These  random  numbers  represent  the  fluctuations  in  the  received 
radar  signal  that  are  caused  by  the  interference,  at  the  antenna  port,  of  the  waves  backscattered 
from  the  different  scattering  centres  of  the  target.  The  statistical  model  can  generate  random 
numbers  appropriate  for  targets  that  have  the  properties  of  the  four  Swerling  models  (see 
Skolnik  [5],  p.  2.21).  The  Swerling  target  type  is  specified  as  an  input  parameter  of  the  target 
simulation  (see  Table  3). 

6.1  Statistical  fluctuations  for  Swerling  1  Targets 

The  Swerling  1  model  describes  a  target  whose  received  signal  voltage  fluctuates  slowly  from 
scan  to  scan  according  to  a  Rayleigh  probability  density  function.  This  model  can  be  used  to 
describe  a  target  with  many  scattering  centres  of  equal  radar  cross  section.  The  relative 
positions  of  these  scattering  centres  with  respect  to  the  radar  do  not  change  significantly  over 
the  time  period  of  the  pulse  repetition  interval,  but  do  change  significantly  over  the  time 
period  between  consecutive  dwells  (i.e.,  the  scan  time). 

A  random  number,  representing  a  statistical  fluctuation,  is  calculated  for  all  of  the  I  and  Q 
values  generated  for  the  entire  scenario.  The  individual  signal  samples  are  represented  as 
\afjn  and  Q„^  where  a  =  0, 1,2,...,  NbNd  -1  is  an  index  that  represents  the  transmitted  burst 

of  the  scenario,  ft  is  an  index  that  represents  the  range  bins  that  are  simultaneously  affected  by 
the  received  signals  corresponding  to  a  particular  burst,  and  n  =  0, 1, 2, ... ,  Nca  - 1  is  an  index 
representing  the  pulses  in  a  burst. 

In  the  Swerling  1  case  four  zero-mean  Gaussian  random  numbers  are  generated  for  each 
complex  signal  sample  as 


f  >  _ 

“*'V2 
f'  -A. 

m 


(46) 
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where  the  xa,  x'a,  x”a,  x"  are  Gaussian  distributed  random  numbers  having  a  mean  of  zero 
and  a  standard  deviation  of  one.  Fluctuations  corresponding  directly  to  the  I  and  Q  values  are 
then  calculated  as 

=[(f«/J2  +fo/J2]  "cos 

(47) 

This  approach  ensures  that  the  fluctuations  corresponding  to  the  quadrature  samples  are  90° 
out  of  phase  with  the  fluctuations  corresponding  to  the  in-phase  samples. 

Equation  (46)  specifies  that  for  a  given  a,  fa ^  is  the  same  for  all  /?  and  n.  That  is,  for  any 

given  burst,  the  value  of  the  random  number  representing  the  fluctuations  is  the  same  for  all  I 
and  Q  values  in  all  range  bins  affected  by  the  returns  corresponding  to  that  particular  burst. 
The  fluctuation  for  the  n*  signal  sample  in  each  range  bin  affected  by  a  particular  burst  must 
be  the  same  since  the  target  motion  will  be  insignificant  over  the  time  period  between  range 
samples,  i.e.,  the  time  required  for  the  radar  signal  to  travel  the  distance  of  the  range 
resolution. 

The  fluctuations  calculated  using  Equation  (47)  are  appropriate  for  a  Swerling  1  target  that  is 
viewed  by  a  multiple  burst  waveform  in  which  the  transmitted  frequencies  change  from  burst 
to  burst.  In  this  case  the  signal  voltage  received  by  the  target  will  fluctuate  faster  than  the  rate 
that  is  defined  for  a  Swerling  1  target.  This  is  because  the  interference  of  the  waves 
backscattered  from  the  scattering  centres  will  vary  as  a  function  of  the  transmitted  frequency. 

In  general,  the  target  simulation  does  not  require  the  frequency  to  change  from  burst  to  burst. 
Thus,  for  a  Swerling  1  target,  the  random  fluctuations  that  are  calculated  in  Equation  (47) 
must  be  adjusted  if  the  frequency  remains  constant  for  any  two  consecutive  bursts  of  a  dwell. 
Consequently,  for  each  dwell  the  frequency  of  burst  0  is  first  compared  with  the  frequency  of 
burst  1.  If  they  are  the  same  then  the  random  fluctuations,  f^  and  f^  ,  corresponding  to 
burst  1  are  set  equal  to  those  corresponding  to  burst  0.  Next  the  frequencies  of  burst  1  and  2 
are  compared.  If  they  are  equal,  the  f^  and  f^  for  burst  2  are  set  equal  to  those  for 
burst  1 .  This  process  is  continued  until,  finally,  the  frequencies  of  burst  ( Nb  -  2)  and 
( Nb  -l)  are  compared.  If  they  are  equal,  the  f^  and  f^  for  burst  (Nb  -l)  are  set  equal  to 
those  for  burst  { Nb  -  2). 

Fluctuations  corresponding  to  the  power  spectral  coefficients  are  calculated  as 

f  Sm  =  fe,  )2  +  te  )2  i  =  n  =  0, 1, 2, . . .,  Nc  - 1  (48) 
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where  the  index  n,  used  in  the  time  domain,  has  been  replaced  with  i,  which  is  used  in  the 
frequency  domain  to  represent  the  coefficients  of  the  power  and  phase  spectra.  The  are 

exponentially  distributed  random  numbers  having  a  mean  of  one.  These  values  are  then 
filtered  with  the  output  of  the  missile  model  (see  Figure  2  and  Equation  (4)).  A  mean  value  of 
one  ensures  that  the  filtering  process  does  not  alter  the  power  of  the  signal  that  is  output  by 
the  missile  model. 

As  stated  in  Section  2,  the  phase  spectrum  of  the  final  signal  for  a  Swerling  1  target  is  taken 
directly  from  the  output  of  the  missile  model.  Hence,  the  statistical  model  does  not  calculate 
phase  fluctuations  in  this  case. 

6.2  Statistical  fluctuations  for  Swerling  2  Targets 

The  Swerling  2  model  describes  a  target  whose  received  signal  voltage  fluctuates  rapidly  from 
pulse  to  pulse  according  to  a  Rayleigh  probability  density  function.  As  in  the  case  of  the 
Swerling  1  target,  this  model  is  used  to  describe  a  target  with  many  scattering  centres  of  equal 
radar  cross  section.  However,  in  this  case  the  relative  positions  of  the  scattering  centres  with 
respect  to  the  radar  do  change  significantly  over  the  time  period  of  the  pulse  repetition 
interval. 

For  a  Swerling  2  target  the  statistical  model  generates  four  random  numbers  for  each  complex 
signal  sample,  i.e., 


where  the  xm ,  x'm ,  ,  x*  are  Gaussian  distributed  random  numbers  having  a  mean  of 

zero  and  a  standard  deviation  of  one.  The  fluctuations  corresponding  to  the  I  and  Q  signal 
samples  are  then  calculated  using  Equation  (47).  In  this  case  the  resulting  ila^  and  f^,  are 

zero-mean  Gaussian  random  numbers  having  a  standard  deviation  of  ^Nca/2  . 

Equation  (49)  specifies  that  unique  random  numbers  are  generated  for  each  signal  sample  of 
each  burst.  However,  as  in  the  case  of  a  Swerling  1  target,  the  random  numbers  are  the  same, 
for  any  particular  signal  sample,  for  each  of  the  multiple  range  bins  that  are  simultaneously 
affected  by  the  returns  from  any  one  burst. 
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Fluctuations  corresponding  to  the  coefficients  of  the  power  and  phase  spectra  are  calculated 
by  DFT,  i.e.,  for  each  a  (burst)  and  for  each  /?  (range  bin)  the  fluctuations  are  calculated  as 


where 


fMRddF  +  Wd)2  f*  =  tan"1 


( Imj 

rfxi\ 

if< 

rfxl 
14;  J) 

(50) 


fx  =(f0x,f1x,f2x,...,fJ_1)=DFT{fI  +jfQ} 

f1  =(f0I,f1I,f (51) 

fQ=(foQ.fiQ.faQ.-.f^) 


The  power  spectral  fluctuations,  f^{ ,  will  be  exponentially  distributed  random  numbers 
having  a  mean  of  one  if  the  DFT  is  defined  with  the  following  normalization 


f(«): =  IT  £fl exp(_  j2/vu  yl 

™ c  y=0 

(52) 

)exp(-  j2nuy/ Nc) 

u=0 


The  power  spectral  fluctuations  are  used  to  filter  the  output  of  the  missile  model  (see 
Equation  (4)).  The  phase  fluctuations,  f^( ,  will  be  random  numbers  uniformly  distributed 

between  -n  and  n.  For  a  Swerling  2  target  these  numbers  are  used  directly  as  the  coefficients 
of  the  phase  spectra. 

6.3  Statistical  fluctuations  for  Swerling  3  Targets 

The  Swerling  3  model  describes  a  target  whose  received  signal  voltage  fluctuates  slowly  from 
scan  to  scan  according  to  a  %2  distribution  of  degree  4  (see  Nathanson  [10],  p.  167).  This 
model  can  be  used  to  describe  a  target  that  consists  of  one  dominant  scattering  centre  and 
many  other  scattering  centres  of  smaller  radar  cross  section.  The  relative  positions  of  these 
scattering  centres  with  respect  to  the  radar  do  not  change  significantly  over  the  time  period  of 
the  pulse  repetition  interval,  but  do  change  significantly  over  the  time  period  between 
consecutive  dwells. 

For  a  Swerling  3  target  the  statistical  model  generates  the  individual  fluctuations 
corresponding  to  the  power  spectral  coefficients  as 
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f|s/  =^[(;c«)2  +  ( x'a )2  +  (XaT  +  U*)2] 


(53) 


where  the  xa ,  x'a,  x",  x”a  are  again  Gaussian  distributed  random  numbers  having  a  mean  of 
zero  and  a  standard  deviation  of  one.  The  resulting,  f^, ,  are  degree  4  x2  distributed  random 
numbers  having  a  mean  of  one. 

Equation  (53)  specifies  that  a  unique  random  number  is  generated  for  each  burst.  As  in  the 
case  of  a  Swerling  1  target,  these  numbers  must  be  adjusted  if  the  transmitted  frequency  is  the 
same  for  any  two  consecutive  bursts  of  a  dwell.  Once  the  f^j,  have  been  appropriately 

adjusted  (see  Section  6.1),  they  are  used  to  filter  the  spectral  power  coefficients  that  are  output 
by  the  missile  model  (see  Equation  (4)). 

As  in  the  case  of  a  Swerling  1  target,  no  fluctuations  are  generated  for  the  coefficients  of  the 
phase  spectra  corresponding  to  a  Swerling  3  target. 

6.4  Statistical  fluctuations  for  Swerling  4  Targets 

The  Swerling  4  model  describes  a  target  whose  received  signal  voltage  fluctuates  rapidly  from 
pulse  to  pulse  according  to  a  x2  distribution  of  degree  4.  As  in  the  case  of  the  Swerling  3 
target,  this  model  can  be  used  to  describe  a  target  that  consists  of  one  dominant  scattering 
centre  and  many  other  scattering  centres  of  smaller  radar  cross  section.  However,  in  this  case 
the  relative  positions  of  the  scattering  centres  with  respect  to  the  radar  do  change  significantly 
over  the  time  period  of  the  pulse  repetition  interval. 

The  individual  fluctuations,  f^, ,  corresponding  to  the  power  spectral  coefficients 
corresponding  to  a  Swerling  4  target  are  calculated  as 

f  Ipi  =^[(*«/)2  +  (*«/  )2  +  {x'ai  Y  +  ( C  Y  ]  (54) 

In  this  case,  unique  random  numbers  are  generated  for  each  spectral  coefficient  of  each  burst. 
These  x2  distributed  numbers  (mean  of  one)  are  then  used  to  filter  the  output  of  the  missile 
model  (see  Equation  (4)). 

The  fluctuations  that  will  represent  the  coefficients  of  the  phase  spectra  are  calculated  as 

fJ,=2*UaI-jr  (55) 

where  the  Uai  are  random  numbers  uniformly  distributed  between  zero  and  one. 
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7  Example  Output 


This  section  presents  example  data  generated  by  the  target  simulation.  The  data  are  presented 
in  the  form  of  comparisons  that  are  intended  to  demonstrate  the  simulation’s  ability  to 
generate  realistic  target  signals.  Specifically,  example  output  generated  by  the  target 
simulation  are  compared  in  terms  of  propagation  factor  and  received  power  with  the  outputs 
of  TERPEM  (TERrain  Parabolic  Equation  Model)  and  ADAPT_MFR  (ADAPTive 
MultiFuntion  Radar  Simulation),  respectively. 

“TERPEM  is  a  propagation  modelling  tool  for  assessing  the  effects  of  atmospheric  refraction 
and  terrain  diffraction  on  radar  and  communications  systems  in  the  frequency  range  30  MHz 
to  100  GHz”  [11].  Since  the  time  of  writing  of  reference  [1 1]  the  lower  frequency  limit  of 
TERPEM  has  been  extended  to  1  MHz  (Levy  and  Craig  [12]).  As  part  of  the  work  to  extend 
TERPEM’s  capabilities,  the  TERPEM  output  was  evaluated  against  the  output  generated  by 
GRWAVE  (Ground-WAVE  propagation  program,  Gill  [13]).  Excellent  agreement  between 
the  outputs  of  the  two  codes  was  achieved  (Levy  and  Craig  [12]).  The  calculations  performed 
by  GRWAVE  are  based  upon  the  theory  of  Rotheram,  [14, 15],  which  represents  a  completely 
different  approach  to  propagation  calculations  than  that  used  by  TERPEM.  The  TERPEM 
output  has  also  been  compared  with  the  output  of  MLAYER.  Again,  excellent  agreement  was 
achieved  [M.  Levy,  private  communication].  MLAYER  was  developed  by  the  Naval  Ocean 
Systems  Centre  for  calculating  the  signal  levels  of  electromagnetic  waves  propagating  in  a 
multilayer  tropospheric  waveguide  environment  over  sea  water  (Yeoh  [16]).  This  software  is 
based  upon  a  full  mode  solution  to  Maxwell’s  equations  and  is  often  used  as  a  baseline  for 
testing  propagation  calculations.  Good  agreement  has  also  been  achieved  (Thomson  [17]) 
between  the  output  of  TERPEM  and  the  outputs  of  AREPS  (Advanced  Refractive  Effects 
Prediction  System,  Patterson  [18])  and  NEMESiS  (Naval  ElectroMagnetic  Environment 
Simulation  System).  Both  of  these  codes  use  parabolic  equation  methods  similar  to  those 
used  in  TERPEM.  NEMESiS  is  based  upon  the  approach  presented  by  Barrios  [19].  Thus, 
TERPEM  is  considered  to  be  accurate  for  the  conditions  of  interest  and  to  represent  a  good 
standard  for  evaluation  of  the  target  simulation  in  terms  of  propagation  factor. 


The  pattern  propagation  factor,  F,  is  defined  “as  the  ratio  of  the  field  strength  that  is  actually 
present  at  a  point  in  space  to  that  which  would  be  present  if  free-space  propagation  had 
occurred  with  the  antenna  beam  directed  toward  the  point  in  question”  (Long  [20],  pill).  In 
terms  of  the  notation  used  in  this  report,  the  TERPEM  propagation  factor  calculations 
consider  the  radar  equation  to  take  the  following  form 


P,G„,GjG[^a, 


(4») 


3  4 

r 


b  p4 

rT 


(56) 


where  Pr  is  the  power  received  from  a  single  radar  pulse  that  is  backscattered  by  a  target  at 
range  r.  By  this  definition  Fj  includes  the  angular  dependence  of  the  antenna  gain  function. 
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The  propagation  factor  is  calculated  from  the  output  I  and  Q  data  generated  by  the  target 
simulation  by  first  determining  the  received  power  as  a  function  of  range.  An  average 
received  power  for  a  particular  range  bin  is  calculated  as 


(57) 


A  propagation  factor  value,  F4 ,  corresponding  to  this  range  bin  is  then  calculated  by 
substituting  Equation  (57)  into  Equation  (56),  i.e., 


F'(r,)  = 


(4ff) 


3  P.  r„4 


P,  Gsys  Gf  G 
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£<J, 


(58) 


These  calculations  have  been  performed  for  the  scenario  defined  in  Tables  1, 2,  and  3.  The 
resulting  propagation  factor  values  (scaled  by  Pc )  are  shown  as  the  green  curve  in  Figure  12. 

However,  it  is  expected  that  these  F4  values  will  not  agree  well  with  the  Fj  values 
calculated  by  TERPEM  (blue  curve  of  Figure  12)  because  of  differences  in  the  form  of  the 
radar  equation  used  by  TERPEM  and  the  target  simulation. 

The  free  space  (path  1)  radar  equation  used  by  the  target  simulation  can  be  written  as 


Propagation  at  20  m  Above  Sea  Level  for  Burst  0 


Figure  12.  Comparison  of  the  propagation  factors  calculated  by  the  target  simulation  and  TERPEM. 
The  modelled  scenario  is  represented  by  the  parameters  listed  in  Tables  1, 2,  and  3.  The  Ts,  f„  and  Nc 
parameters  for  this  case  represent  burst  Oof  a  four  burst  dwell  that  is  repeated  continuously  throughout 

the  flight  time  of  the  missile.  The  blue  curve  shows  Fr4  vs.  range.  The  red  curve  shows  F'4  vs. 
range.  The  green  curve  shows  F4  / Pc  vs.  range. 
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where  G2N(0t,  <f>tAc, &)  =  G2/{G2f  G>)  is  the  normalized  antenna  gain  pattern.  This 

equation  determines  the  power,  corresponding  to  a  range  bin  at  rc ,  that  is  received  from  a 

single  transmitted  pulse.  The  presence  of  the  earth  can  be  incorporated  into  this  radar 
equation  by  adding  a  pattern  propagation  factor,  which  gives 


„  /-  v  PA  *2  Gsys  G)  G2  W2(rt,  rc)ab 
— (4 * 


(60) 


Thus,  propagation  factor  values  calculated  using  Equations  (57)  and  (58)  would  incorporate 
the  effects  of  pulse  compression,  attenuation  by  precipitation,  and  range  weighting,  which  are 
not  considered  by  TERPEM.  Therefore,  to  achieve  a  better  comparison,  the  scenario  defined 
by  Tables  1,  2,  and  3  is  executed  by  the  target  simulation  again,  with  the  attenuation  and  range 
weighting  calculations  switched  off.  In  addition,  since  TERPEM  only  performs  calculations 
in  a  two  dimensional  plane  (rangd  and  height),  the  input  parameter  of  the  target  simulation 
describing  the  target  azimuth  is  set  equal  to  the  azimuth  of  the  radar  beam  axis 
(<ps  =  0c  =  130° ).  TERPEM  also  does  not  consider  statistical  fluctuations.  Hence,  the 
Swerling  number  input  parameter  is  set  to  zero,  which  corresponds  to  no  statistical 
fluctuations.  Using  these  considerations,  average  received  powers,  P' ,  are  calculated  from 
the  I'  and  Q'  outputs  as 


f;fc)=TVSl(r-)5+(Q;)!l  (61) 

A  propagation  factor  value  ,  F'4  ,  that  is  comparable  with  that  produced  by  TERPEM  is  then 
calculated  for  each  range  bin  as 


F'*(rc)  = 


(4;r)3  P'  r4 


P'G^GjG 


2  12 


(62) 


The  red  curve  plotted  in  Figure  12  shows  the  range  dependence  of  these  F'4  values. 
Agreement  to  within  3  dB  of  the  Fr4  values  produced  by  TERPEM  (blue  curve)  is  now 
achieved  at  all  range  bins,  even  in  the  multipath  nulls. 

The  waveform  parameters,  Ts,f,  and  Nc  corresponding  to  the  case  of  Figure  12  (see  Table  2) 
describe  one  burst  of  a  four  burst  dwell  waveform.  The  target  simulation  generated  the 
plotted  data  by  transmitting  this  dwell  waveform  repeatedly  throughout  the  flight  time  of  the 
missile.  The  waveform  parameters  describing  all  four  bursts  of  the  dwell  are  shown  in 
Table  4.  Similar  comparison  plots  between  the  propagation  factors  calculated  by  TERPEM 
and  those  derived  from  the  target  simulation  for  bursts  1,  2,  and  3  are  shown  in  Figures  13, 14, 
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Table  4.  Waveform  parameters. 


BURST  NUMBER 

Ts(ps) 

/,(GHz) 

Nc 

0 

240.7 

10.93 

26 

1 

221.5 

10.65 

29 

2 

212.8 

9.74 

33 

3 

206.4 

8.955 

37 

and  15.  Similarly  good  agreement  between  the  target  simulation  and  TERPEM  is  achieved 
for  each  burst. 


As  stated  in  Section  2,  the  target  simulation  will  not  correctly  calculate  the  propagation  of  the 
radar  signal  into  the  “intermediate  region”.  Blake  [4]  (p.  279)  gives  quantitative  expressions, 
based  on  grazing  angle,  that  approximate  the  region  of  validity  of  the  geometrical  optics 
approach  used  in  the  target  simulation.  However,  these  expressions  are  not  exact.  A  better 
method  of  testing  whether  or  not  the  geometrical  optics  approach  is  valid  may  be  to  compare 
the  calculated  results  with  those  of  TERPEM  as  has  been  done  in  Figures  12, 13, 14,  and  15. 
In  these  cases  the  target  simulation  agrees  well  with  TERPEM  and  can  be  considered 
accurate. 

Figure  16  shows  an  example  for  which  the  target  simulation  becomes  inaccurate.  In  this  case, 
the  parameters  corresponding  to  Figure  15  have  been  used  except  for  the  antenna  height, 
which  has  been  lowered  to  5  m.  Lowering  the  antenna  height  causes  the  radar  horizon  to 
decrease.  Thus,  the  target  simulation  calculations  extend  into  the  intermediate  region.  The 
propagation  factor  calculated  by  the  target  simulation  now  deviates  from  that  calculated  by 
TERPEM  at  far  ranges.  Therefore,  the  target  simulation  calculations  should  only  be  used,  for 


Propagation  at  20  m  Above  Sea  Level  for  Burst  1 


Range  (km) 

Figure  13.  Comparison  of  the  propagation  factors  calculated  by  the  target  simulation  and  TERPEM. 
The  target  simulation  data  corresponds  to  the  burst  1  waveform  parameters.  The  blue  curve  shows 

Fj  vs.  range.  The  red  curve  shows  F'4  vs.  range.  The  green  curve  shows'  F4  /  Pc  vs.  range. 
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Propagation  at  20  m  Above  Sea  Level  for  Burst  2 


Range  (km) 

Figure  14.  Comparison  of  the  propagation  factors  calculated  by  the  target  simulation  and  TERPEM. 

The  target  simulation  data  corresponds  to  the  burst  2  waveform  parameters.  The  blue  curve  shows 

F?  vs.  range.  The  red  curve  shows  FH  vs.  range.  The  green  curve  shows  F4 / Pc  vs.  range. 

the  case  of  Figure  16,  at  ranges  less  than  ~  22.5  km.  This  method  of  assessing  the  accuracy  of 
the  geometrical  optics  approach  used  in  the  target  simulation,  i.e.,  comparison  with  TERPEM, 
leads  to  the  obvious  recommendation  that  the  target  simulation  should  be  upgraded  to  use 
TERPEM  for  propagation  calculations. 

ADAPT_MFR  is  the  latest  in  a  line  of  comprehensive  computer  programs  (developed  under 
contract  for  the  Defence  R&D  Canada  -  Ottawa)  designed  to  evaluate  the  detection  capability 
of  multifunction  radars  when  operating  against  anti-ship  missiles  (Gauthier  et  al.  [21], 
Gauthier  [22]).  The  simulated  scenarios  can  include  up  to  five  missiles  (sea  skimmers  or  high 
divers)  and  up  to  five  collocated  multifunction  radars  operating  simultaneously  in  littoral 
regions  in  the  presence  of  jamming.  These  radars  are  optimized  for  horizon  search  but  can 
perform  any  user-defined  search  pattern,  including  cueing,  confirmation  and  tracking.  As  a 
consequence  of  the  larger  scope  of  the  modelled  scenarios,  relative  to  the  target  simulation, 
the  ADAPT_MFR  calculations  are  performed  with  lower  fidelity,  e.g.,  the  simulation  does  not 
calculate  the  complex  components  of  the  received  signal  voltage.  It  does,  however,  calculate 
the  received  target  powers  as  a  function  of  range  using  TERPEM  output  for  the  propagation 
factor.  Since  these  calculations  are  performed  in  a  different  manner  than  those  of  the  target 
simulation,  a  favourable  comparison  between  the  two  simulation  outputs  should  increase 
confidence  in  the  accuracy  of  the  outputs  produced  by  both  simulators. 

The  burst  0  target  power  as  calculated  by  the  target  simulation  using  the  parameters  of 
Tables  1,  2,  and  3  is  plotted  as  a  function  of  range  (green  curve)  in  Figure  17.  The  target 
power  as  calculated  by  ADAPT_MFR  using  equivalent  parameters  is  also  plotted  (blue  curve) 
in  Figure  17.  As  was  the  case  for  the  propagation  factor  comparisons  of  Figures  12, 13, 14, 
and  15,  the  main  difference  between  these  two  curves  can  be  accounted  for  by  factors  that  are 
not  considered  by  both  simulations.  In  this  case,  ADAPT_MFR  only  calculates  the  mean 
target  power.  Statistical  fluctuations  of  the  received  target  power  are  only  considered  in  the 
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Propagation  at  20  m  Above  Sea  Level  for  Burst  3 


Range  (km) 

Figure  15.  Comparison  of  the  propagation  factors  calculated  by  the  target  simulation  and  TERPEM. 

The  target  simulation  data  corresponds  to  the  burst  3  waveform  parameters.  The  blue  curve  shows 

vs.  range.  The  red  curve  shows  F'4  vs.  range.  The  green  curve  shows  F4 / Pc  vs.  range. 

probability  of  detection  calculations.  In  addition,  ADAPT_MFR  does  not  explicitly  consider 
the  range  weighting  applied  to  targets  at  different  locations  within  the  measurement  volume. 
Thus,  the  output  of  the  target  simulation  is  regenerated  with  the  statistical  fluctuations  and 
range  weighting  calculations  turned  off.  These  adjusted  data,  plotted  as  the  red  curve  in 
Figure  17,  agree  very  well  with  the  results  calculated  by  ADAPTJMFR.  The  small 
differences  between  the  adjusted  data  produced  by  the  target  simulation  and  the  data  produced 
by  ADAPT_MFR  for  equivalent  conditions  are  almost  entirely  due  to  the  differences  found 
between  the  propagation  factors  calculated  by  TERPEM  and  the  target  simulation.  Similar 
plots  corresponding  to  bursts  1,  2,  and  3  are  shown  in  Figures  18, 19,  and  20.  Similarly 
favourable  comparisons  are  achieved  in  each  case. 
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Propagation  at  5  m  Above  Sea  Level  for  Burst  3 


Figure  16.  Comparison  of  the  propagation  factors  calculated  by  the  target  simulation  and  TERPEM. 
The  parameters  of  the  modelled  scenario  are  the  same  as  that  corresponding  to  Figure  15,  except  for 

antenna  height,  which  was  set  to  5  m.  The  blue  curve  shows  F4  vs.  range.  The  red  curve  shows  F'4 
vs.  range.  The  green  curve  shows  F4 / Pc  vs.  range. 


Burst  0  Target  Power 


Figure  17.  Comparison  of  the  burst  0  target  powers  generated  by  the  target  simulation  and 
ADAPT_MFR  for  the  scenario  defined  in  Tables  1, 2,  and  3.  The  green  cun/e  shows  the  standard 
output  of  the  target  simulation.  The  red  curve  shows  the  data  generated  by  the  target  simulation  without 
range  weighting  or  statistical  fluctuations.  The  blue  curve  shows  the  equivalent  target  powers  generated 

by  ADPT_MFR. 
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Power  (dBm) 


Burst  1  Target  Power 


10  15  20  25  30 

Range  (km) 


Figure  18.  Comparison  of  the  burst  1  target  powers  generated  by  the  target  simulation  and 
ADAPT_MFR  for  the  scenario  defined  in  Tables  1,  2,  3,  and  4.  The  green  curve  shows  the  standard 
output  of  the  target  simulation.  The  red  curve  shows  the  data  generated  by  the  target  simulation  without 
range  weighting  or  statistical  fluctuations.  The  blue  curve  shows  the  equivalent  target  powers  generated 

by  ADPT_MFR. 
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Power  (dBm) 


Figure  19.  Comparison  of  the  burst  2  target  powers  generated  by  the  target  simulation  and 
ADAPT_MFR  for  the  scenario  defined  in  Tables  1,  2,  3,  and  4.  The  green  curve  shows  the  standard 
output  of  the  target  simulation.  The  red  curve  shows  the  data  generated  by  the  target  simulation  without 
range  weighting  or  statistical  fluctuations.  The  blue  curve  shows  the  equivalent  target  powers  generated 

by  ADPT_MFR. 
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Burst  3  Target  Power 


Figure  20.  Comparison  of  the  burst  3  target  powers  generated  by  the  target  simulation  and 
ADAPT_MFR  for  the  scenario  defined  in  Tables  1,  2,  3,  and  4.  The  green  curve  shows  the  standard 
output  of  the  target  simulation.  The  red  curve  shows  the  data  generated  by  the  target  simulation  without 
range  weighting  or  statistical  fluctuations.  The  blue  curve  shows  the  equivalent  target  powers  generated 

by  ADPT_MFR. 
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8  Summary  and  Conclusions 


This  report  describes  a  high  fidelity  simulation  of  the  complex  voltage  signals  produced  by 
the  receiver  of  a  shipbome  horizon  scanning  phased  array  radar  when  observing  an 
approaching  sea-skimming  anti-ship  missile  in  adverse  weather  conditions.  The  output  of  this 
simulation  is  to  be  used  together  with  the  output  of  high  fidelity  clutter  and  radar  signal 
processor  simulations  to  investigate  radar  detection  performance  (e.g.,  see  Thomson  [3]). 
Weather  effects  are  only  considered  in  terms  of  attenuation  of  the  signals  received  from  the 
target.  The  attenuation  models  used  limit  the  applicability  of  this  simulation  to  radar 
frequencies  in  the  X,  C,  or  S-band. 

The  modelled  target  moves  directly  toward  the  observing  radar  at  constant  speed  along  a 
flight  path  of  constant  altitude  and  azimuth.  The  received  target  signals  can  incorporate 
statistical  fluctuations  corresponding  to  the  any  of  the  four  Swerling  target  types. 

The  environment  is  modelled  as  a  stratiform  precipitation  field  above  a  rough  spherical  sea 
surface,  exactly  as  was  done  in  Thomson  [1],  Refraction  of  the  radar  signals  by  the 
atmosphere  is  considered  by  using  an  effective  earth  radius  approach.  Multipath  interference 
caused  by  the  reflection  of  the  radar  signal  by  the  sea  surface  is  calculated  geometrically. 
Diffraction  effects  are  ignored.  Thus,  the  applicability  of  this  target  simulation  is  limited  to 
the  “interference  region”  (see  Blake  [4],  p  279)  of  atmospheres  having  an  index  of  refraction 
that  changes  linearly  with  height. 

The  modelled  radar  system  transmits  a  multiple  burst  waveform  in  which  the  frequency,  pulse 
repetition  time,  and  number  of  pulses  per  burst  can  all  change  from  burst  to  burst  and  dwell  to 
dwell.  The  motion  of  the  target  from  pulse  to  pulse  is  considered.  Range  weighting  of  the 
target  within  the  resolution  volume  is  also  considered.  The  specific  forms  of  the  range 
weighting  and  antenna  gain  functions  used  in  the  target  simulation  represent  upgrades  to  the 
models  used  in  the  weather  clutter  simulation  as  described  by  Thomson  [1],  These  upgrades 
have  now  also  been  included  in  the  weather  clutter  simulation. 

The  target  simulation  was  evaluated  by  comparing  its  output  in  terms  of  pattern  propagation 
factor  and  received  target  power  with  the  outputs  of  TERPEM  [11]  and  ADAPT_MFR 
(Gauthier  et  al.  [21];  Gauthier  [22],  respectively.  Good  agreement  was  achieved  in  each  case, 
except  where  the  calculations  were  extended  beyond  the  interference  region  into  the 
intermediate  region.  To  increase  the  applicability  of  this  simulation  it  is  recommended  that  it 
be  integrated  with  TERPEM  so  that  propagation  can  be  correctly  considered  at  any  range  in 
any  atmospheric  structure.  Nevertheless,  the  favourable  comparisons  in  the  interference 
region  demonstrate  that  this  simulation  is  accurate,  where  applicable,  and  that  it  can  be  used 
with  confidence  for  its  stated  purpose  of  assessing  radar  detection  performance.  In  addition, 
the  good  agreement  between  the  target  simulation  and  ADAPT_MFR  allows  further 
comparisons  to  be  made  between  the  two  simulations,  e.g.,  in  terms  of  probability  of 
detection,  with  the  knowledge  that  potential  differences  are  not  due  to  modelling  of  the  power 
received  from  the  target. 
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Annex  A:  Reflectance  of  the  Sea  Surface 


The  reflectance  describes  the  ratio  of  reflected  to  incident  energy.  This  quantity  must  be 
calculated  for  the  case  of  a  rough  spherical  sea  surface.  Typically,  the  roughness  of  the  sea  is 
considered  to  be  caused  by  wind  and  swell,  and  is  characterized  in  terms  of  a  sea  state. 
However,  precipitation  impacts  at  the  reflection  point  may  also  have  an  effect  on  the 
reflectance.  There  are  no  works  known  to  the  author  that  consider  the  effect  of  a  precipitation 
disturbed  sea  surface  on  the  reflected  wave.  Backscatter  from  a  rain  disturbed  sea  surface  has, 
however,  been  investigated.  In  these  studies  the  time  history  of  the  splash  structure  resulting 
from  a  raindrop  impact  has  been  investigated  (Hansen  1986  [23]).  Three  main  scattering 
structures  have  been  found.  These  are  the  “crown”  and  “stalk”  of  the  splash  as  well  as  the 
“ring  waves”  that  propagate  outward  along  the  sea  surface  from  the  impact  location.  Since 
the  crown  and  stalk  features  are  above  the  sea  surface,  the  resulting  backscatter  may  be  more 
similar  to  that  of  rain  backscatter  than  that  of  sea  backscatter.  Therefore,  in  this  study  these 
features  are  considered  to  result  in  a  negligible  increase  in  rain  rate  at  the  sea  surface,  which 
will  cause  a  negligible  increase  in  attenuation  of  the  forward-scattered  wave.  It  is  also 
assumed  that  the  small-amplitude  ring  waves  do  not  introduce  significant  additional  structure 
to  the  already  existing  structure  caused  by  gravity  waves  and  capillary  waves.  Thus,  it  is 
assumed  that  disturbance  of  the  sea  surface  by  impacting  precipitation  does  not  significantly 
effect  the  reflectance. 

Rain  will  also  have  a  significant  effect  on  the  temperature  and  salinity  of  the  sea  surface, 
depending  on  the  rain  rate  and  duration  of  the  rainfall  (Katsaros  and  Buettner  1969  [24]). 
However,  the  reflectance  of  the  sea  surface  is  not  very  sensitive  to  changes  in  temperature  and 
salinity.  Therefore,  this  effect  is  also  neglected  in  the  reflectance  calculations. 

The  reflectance,  p,  is  given  by  the  magnitude  of  the  complex  reflection  coefficient,  i.e., 

„  H 

w 

The  reflection  coefficients  for  a  rough  spherical  sea  surface  are  given  by  Skolnik  [5]  (p.  2.38) 
as 


vertical  polarization 
horizontal  polarization 


(63) 


Rv=RwFrFD  vertical  polarization 

(64) 

Rh  =  Rho  Fr  Fd  horizontal  polarization 

where  the  Rm  and  Rh„  are  the  magnitudes  of  the  smooth-flat-earth  reflection  coefficients  for 
vertical  and  horizontal  polarizations,  respectively.  The  specular-reflection  roughness  factor  is 
given  by  Fr ,  and  FD  denotes  the  divergence  factor. 

Skolnik  [5]  (p.  2.40)  gives  the  specular-reflection  roughness  factor  as 

Fr=e~?  J0(j£)  (65) 


44 


DRDC  Ottawa  TR  2002-145 


where  J0  is  the  order  zero  Bessel  function  of  the  first  kind  and 


#  =  2 


2n  Hsin^j^ 

A 


(66) 


The  standard  deviation  of  wave  height,  H,  is  parameterized  in  terms  of  a  sea  state  parameter 
(see  Table  5),  which  is  a  required  input  of  the  simulation. 

The  grazing  angle,  y/g3 ,  is  calculated  using  the  path  3  geometry  depicted  in  Figure  11.  It  has 

been  shown  that  a  cubic  equation  in  the  arc  length,  au  can  be  written  in  terms  of  ha,  Re,  ht,  and 
the  arc  length  a  (see  Freehafer  et  al.  [25],  p.  1 13),  i.e., 

2flj3  -  3 aal  +  [a 2  - 2 Re (ha  +ht )]  aL  +  2 Rehaa  =  0  (67) 


The  solution  to  Eq.  (67)  is  given  by 


<*1 


a 

=  —  +  q,  cos 
2  1 


'<h+n' 
{  3 


(68) 


where 


Using  the  law  of  sines  on  the  outer  triangle  of  Figure  1 1  the  arc  length  a  can  be  written  as 


a  =  R„  sin 


-l 


rt  sin0( 
R„  +  h 


t  J 


(70) 


Once  a\  is  determined,  the  angle  /?3  =  «[  / Re  can  be  used  to  calculate  the  range  to  the 
reflection  point,  lm ,  by  applying  the  law  of  cosines  to  the  left  triangle  of  Figure  11,  i.e., 


Table  5.  Relationship  of  sea  state  to  standard  deviation  of  wave  height. 


SEA  STATE,  Hs 

H(m) 

SEA  STATE,  Hs 

H(m) 

0 

0.0 

4 

0.4953 

1 

0.0762 

5 

0.762 

2 

0.1524 

6 

1.2192 

3 

0.3048 

7 

2.286 
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<,„■  =fc  +  («,  +  K)2 -M.(K  +  *.)««a]^  (71) 

The  range  from  the  reflection  point  to  the  target  is  similarly  calculated  by  applying  the  law  of 
cosines  to  the  right  triangle  of  Figure  1 1,  i.e., 

to,  =fc  +  («,  +h,Y  +*, )cos(A  -A)]l 

where  /?,  =  a/Re .  The  range  £tV  is  used  to  calculate  the  polar  angle  of  the  reflection  point 
by  applying  the  law  of  sines  to  the  left  triangle,  i.e., 


0t3  =  it-  sin 


K sinA^I 

**  J 


(72) 


Finally,  the  grazing  angle  can  be  determined  from  /?3  and  0,3 ,  i.e., 


¥g3  ~  ®t3  03  2 


(73) 


The  divergence  factor  describes  the  weakening  of  the  signal  power  density  caused  by 
reflection  from  a  spherical  surface.  It  is  calculated  as 


f  Ji+Miz£l> 

[  aRe  sini/fgi 


Although  the  atmosphere  was  considered  to  be  an  absorbing  medium  for  the  attenuation 
calculations,  it  will  be  assumed  to  be  non-absorbing  when  calculating  the  reflection 
coefficients.  For  this  case,  the  equations  of  Fowles  [26]  (p.  166-167)  can  be  used  to  write  the 
smooth-flat-earth  reflection  coefficients  as 


-N2sin^3  +ij N2  -  cos2  yrg3 
N2sin^3  +  yN2  -  cos2 

(75) 

sin yrg 3  -  ^N2  -  cos1  iff g3 
sin +  ^N2  -  cos2  iffgi 


where  N  is  the  complex  index  of  refraction  of  sea  water.  For  a  non-magnetic  material, 

N  =  ,  where  er  is  the  complex  dielectric  constant  (or  relative  permittivity).  The  complex 

dielectric  constant  for  sea  water  is  calculated  using  the  equations  of  Klein  and  Swift  [27],  The 
standard  Debye  expression  gives  er,  in  F/m,  as 
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(76) 


j<7 

e  =  e  +  — 2 — - - - - 

\+]2nitrr  2nite0 

where  £L=  is  the  dielectric  constant  at  infinite  frequency,  es  is  the  static  dielectric  constant,  f,  is 
the  transmitted  radar  frequency,  Tr  is  the  relaxation  time,  and  <ris  the  ionic  conductivity. 

Klein  and  Swift  [27]  state  that  choosing  ex  -  4.9  F/m  ±  20%  contributes  negligible  error.  The 
remaining  three  quantities,  es,  rn  and  (Tare  parameterized  in  terms  of  the  sea  surface 
temperature,  Tsea  in  °C,  and  salinity,  s  in  parts  per  thousand,  i.e., 

£s(Tsea’S)  =  £s(Tsea)a£(S,Tsea) 

es (Tsea )  =  87. 134  - 1 .949  x  10-1  Tsea  - 1.276  x10“2t1 +2.491  xlO'4^ 
a£(s,  Tsea)  =  \. 0  + 1.613 x  10‘557;ra  -  3.656 x  10_3s  +  3.21  x  10'V  -  4.232 x  10_V 

*r  (Tsea  ^)  =  tr  (Tsea ,  0)  b£  (s,  Tsea  )  / 

Tr (Tsea ,  0)  =  1 .768  x  1 0_1 1  - 6.086 x  10“137;M  + 1. 104 x  10"14rs2a  -8.111xlO‘17Ts3a  (77) 

be (s,  Tsea )  =  1 .0 +  2.282 x KT5s7;ea  - 7.638 x lO-4*  - 7.76 x HT6 s2  + 1 . 105 x KT8  j3 

<?(Tsea .  s)  =  cr(25,  s)  exp(-  A  Be ) 

A  =  25-7;efl 

Be  =  2.033  xlO~2  +  1.266x10"^  +  2.464xl0'6A2 

-  s(l  .849  x  10'5  -  2.55 1  x  1  O'7  A  +  2.55 1  x  10"8  A2  ) 

<7(25,  s)  =  .s(o.  18252 1-1. 46 192 xlO_3,s  + 2.09324 xlO_5s2  -1. 28205 x10-7s3) 
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